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Abstract 

We consider the nonlinear Klein Gordon Maxwell system derived from the Lagrangian J ( K ~\F IJlv F 111 ' + 
\{{d — ieA) M </>, (d — ieA) M </>) — V(4>) — eA M J^) on four dimensional Minkowski space-time, where is a 
complex scalar field and F^ — 9 M A„ — <9„A M is the electromagnetic field. For appropriate nonlinear 
potentials V, the system admits soliton solutions which are gauge invariant generalizations of the non- 
topological solitons introduced and studied by T.D. Lee and collaborators for pure complex scalar fields. 
In this article we develop a rigorous dynamical perturbation theory for these solitons in the small e limit, 
where e is the electromagnetic coupling constant. The main theorems assert the long time stability of the 
solitons with respect to perturbation by an external electromagnetic field produced by the background 
current J s , and compute their effective dynamics to 0(e). The effective dynamical equation is the 
equation of motion for a relativistic particle acted on by the Lorentz force law familiar from classical 
electrodynamics. The theorems are valid in a scaling regime in which the external electromagnetic fields 
are 0(1), but vary slowly over space-time scales of 0(|), and S = e 1_fc for k £ (0, |) as e -> 0. We work 
entirely in the energy norm, and the approximation is controlled in this norm for times of O(i). 
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1 Statement of results 

1.1 Introduction 

In this article, we are interested in the effective dynamics of a class of solitary wave, or soliton, solutions 
to the nonlinear Klein-Gordon-Maxwell (nl-KGM) equations, in the presence of an external electromagnetic 
field. In this introduction we start by writing down the equations and giving a heuristic statement of, and 
motivation for, our results in $1. 1.31 and t U.1.41 Then, in i jl.2l and M1.31 we provide the necessary background 
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for a precise formulation of the main results - theorems [TO] and [TO] - which appear in i jl.4l These theorems 
are proved in the subsequent sections; a list of notation appears in ill. 1.51 to facilitate reading of the article. 

1.1.1 The equations 

We study the following system of equations, called the nonlinear Klein-Gordon-Maxwell system, or (nl-KGM) 
system, which describe the interaction of a complex scalar field <fr with an electromagnetic field F^ v in the 
presence of an external space-time current JJ B : 

d> i F llv = e{i(t>,^ v (l>)+eJ* 

B^(f> + V'(cj>) = 0. ' 

Here is a complex function on Minkowski space-time R 1+3 , and = <9 M — ieA^ is the covariant derivative 
associated to an electromagnetic potential A^dx^ = Aadt + Ajdx^ with associated field — d^A v — 
d v A^. (The operator D determines an S* 1 connection over R 1+3 whose curvature is —iF.) We use standard 
relativistic notation in which {x^}^Zg are co-ordinates, with greek indices running over {0,1,2,3}, x° = t 
is the time co-ordinate, and {x->}j_± are space co-ordinates with latin indices running over {1,2,3}; the 
Minkowski metric is 

^dx^dx" = dt 2 - (dx 1 ) 2 - (dx 2 ) 2 - (dx 3 ) 2 , 

and is used to raise/lower indices in the usual way. When the spatial part of a space time vector or 1-form 
is considered separately bold face will often be used e.g. x = (x 1 ,^ 2 ,^ 3 ) for clarity. We refer to e as the 
(electromagnetic) coupling constant: for the purposes of this article it is a small positive parameter. The 
current four-vector is of the form 

J B =J B ' v d„ = p B dt+j k B d k 

and is conserved, i.e. 

dtPB + divjs = 0. 

The quantity p B is called the (background) charge density, while js is referred to as the (background spatial) 
current density. Throughout the paper we make the following hypotheses on the nonlinear potential function 
V: 

(HI) Phase invariance: there exists G : R — > M such that V{<j>) = G(\</)\). 

(H2) Positive mass: V{<j>) = +Vi((p) where m > and Vi(0) = -U(\<j>\) is smooth with 17(0) = 

U'(0) = U"(0) = 0. 

(H3) Sub-criticality: the third derivative D^V\ — V 1 satisfies a growth condition |Vj (<fi)\ < c(l + |^| p_3 ), 
for some p G (3, 6). The significance of 6 is that it is the critical Sobolev exponent for the embedding 
H 1 (R i ) L P (R 3 ). 

The function V is subject to a number of additional more specialized hypotheses, which we detail in 
til. 3. 21 in particular to ensure existence and uniqueness of solitons solutions with the properties described in 

El 

1.1.2 Solitons 

The research in this paper is built upon the existence results for solitons in semi-linear wave equations given 
m 0,11] and [12] • These solitons are time-periodic solutions of the nonlinear Klein-Gordon equation 

d fl d»<p + V'(<p) = 0, 

which is obtained by putting e = in ([1]) (i.e. when there is no electromagnetic coupling), and are of the 
form 

0(t,x) = c <w */ w (x). 
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T.D. Lee emphasized that solutions of this type, which he called non-topological solitons, provide a way of 
circumventing the Derrick-Pohozaev non-existence results on static solitons in scalar field theories; see [T5l 
Chapter 7] for a discussion of their properties from the physical point of view. 

It is proved in the references [3], [I] and [22] that, for certain potentials V, solutions of this form exist with 
fuj positive and radial. Also under further conditions these solutions are known to be essentially unique ( |18j ) 
and dynamically stable ( 10, 24 ); see HI. 31 and the appendices for further details. For non-zero values of the 
coupling constant e solutions to {J) of this type have been constructed in [2] [1] directly, using a spherically 
symmetric ansatz, and perturbatively in [161 117) for small e using the e = case as a starting point. For 
small e it is possible to use the information on stability for e = from [24j to prove modulational stability 
of the solitons and their Lorentz boosts, see £11.3.41 and [T7] for details. Much of the same information for 
the e = case will also be used in the present article to study the stability of the solitons when subjected to 
external (background) electromagnetic fields. 

1.1.3 Informal statement of results on interaction of solitons with electromagnetic field 

Our main concern in this article is to understand the interaction of the solitons just described, with an 
external electromagnetic field produced by the space-time current J s . In order to be able prove theorems 
giving precise information on the effect of this field on the soliton, we study ([I]) in a regime determined by 
two small parameters: 

• The electromagnetic coupling constant e = o(l). 

• The external electric and magnetic fields, ~Ei S ext and Bf^, vary over scales which are 0(g), where 
S = o(l). Thus the small parameter 5 is the ratio of the size of the soliton to the length scale over 
which the external field varies. 

The following is an informal version of our main theorems: 

The system ([T]) has solutions which are close, in energy norm, to solitons of the type described 
above and which, in an appropriate scaling regime, move according to the Lorentz force equation: 

j t ( 7 M s u) = eQ s (~E s ext + u x Bf xt ) , (2) 

where the effective mass Mg and charge Qs of the soliton are as in (|60[) and (|6ip . The scaling 
S = e 1 ~ k for k £ (0, ^) ensures that this holds for time intervals of length as e — > 0. 

The precise formulation is in the two theorems stated in ^1.41 

1.1.4 Motivation and related work 

Our interest in this problem stems from the classical, but ongoing, controversy surrounding the classical 
equation of motion for a point charge in an external electromagnetic field. The difficulty arises in attempts 
to account for the "back reaction" of the charge's own field on itself. Attempts to derive an equation of 
motion lead to modifications of the Lorentz force law ([2]) , most notably the Lorentz-Dirac equation ( |21[ 
Equation 9.1]). This equation is third order in time, and is difficult to interpret consistently without some 
further constraint on the type of solution allowed, due to the occurence of runaway solutions and violations 
of causality, (see [6] and [7[ Chapter 28]). Recent discussions of this problem have been given in [9] and the 
books [2U12H]. One natural and well established approach to the problem of making sense of the back reaction 
is to start with a well-posed system of equations in which the point charge is explicitly replaced by a smooth 
bounded charge distribution, the Abraham model, or one of its generalizations like the Lorentz model, for 
example. One can then derive an equation of motion for the charge as an expansion, valid when the size of 
the charge distribution is small (compared with typical length scales set by the external fields), and show 
that this agrees with the Lorentz-Dirac equation at a certain order of approximation - see |14j . In this setting 
it turns out, however, that at the same order the Lorentz-Dirac equation can be approximated by a more 
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conventional equation of motion which seems to be free of interpretational difficulties, (see [2TJ Equation 
9.10], where the name Landau-Lifshitz equation is suggested for this effective equation of motion. The 
Landau-Lifshitz equation, which is second order in time, can be obtained formally from the Lorentz-Dirac 
equation by substituting for the third derivatives the expression obtained by differentiating the ordinary 
Lorentz force law once in time.) 

Our aim in studying soliton motion in the (nl-KGM) system is to attempt a similar analysis using a 
solitonic model for the particle (in place of the Abraham or Lorentz model). Our model has the virtue 
of being, in a very natural way, a Lorentz invariant system which is well posed (and so free of causality 
problems). Unfortunately the calculations required even just to derive the equation of motion for the soliton 
to 0(e) (i.e. the Lorentz force equation @) are long, and further work will be required to calculate additional 
corrections which may be compared with the Lorentz-Dirac equation in appropriate regimes. To achieve this, 
the starting point would be the equation of motion (|116[) for the soliton parameters derived from modulation 
theory. In $4] this equation is computed to highest order (i.e. to 0(e)), and shown to give the Lorentz 
force law. A computation to the next order should give the Landau-Lifshitz equation ([2lJ Equation 9.10]). 
However it seems that some renormalization of the soliton mass and charge will have to be taken into account 
in this computation, and it is possible a refinement of the ansatz ([52)1 will be needed to achieve this. It is to 
be hoped that at least in some simple cases such as one dimensional motion of the soliton in an electric field 
^ixt = (0,0, E(St,Sx)) of fixed direction it will be possible to carry this through, and make a comparison 
with the corresponding specialization of the Landau-Lifshitz equation ([5TJ Equation 9.11]). 

A corresponding theorem to our main result was proved for solitons in interaction with gravitational 
fields in the articles [551 US] • The system treated there (Einstein's equation coupled to a nonlinear klcin- 
Gordon equation) is in many ways more difficult than the one studied here (for example it is quasi- linear). 
Correspondingly, it is possible to carry out a more general analysis for the Klein-Gordon-Maxwell system 
under consideration here: in particular we emphasize that in the present article we are able to work entirely 
with the energy norm throughout (whereas for the Einstein system it was necessary to work with much 
stronger norms). There have also been theorems proved on effective dynamics for solitons moving under a 
potential in the nonlinear Schroedinger equation, see [TJl [5J [S7J . 

1.1.5 Notation 

The following is a list of notations for important objects, with the section in which they are first introduced, 
for reference. 

• L P (R 3 ) is the Lebesgue space of (equivalence classes of) measurable functions with norm ||/||lp = 
Jr3 \ f\ P dx < co, and H k (M. 3 ) is the Sobolev space of (equivalence classes of) measurable functions 
with norm = X)|a|=o ll^"/IU 2 < °°j where d a means the weak partial derivative determined by 
the multi-index a. We say / £ Hh if \f G H k for every smooth, compactly supported \-> an< i 

H 1 = {/ G H\ oc n L 6 : ||V/|| L2 = < oo} . (3) 

Further we define H k to be the intersection of H k and the space of radial functions, i.e. functions of 
1 2 1, and similarly define and H}. 

• Electromagnetic potential A^dx^ = A dt + Ajdx^ , electromagnetic field F^ v = c^A,, — d„A^, and 
covariant derivative B M = <9 M — ieA^: 31.1.11 

• Complex scalar (soliton) field <f> and its self- interaction potential V(4>) — G(\<fi\) subject to hypotheses 
(H1)-(H3): flTXTl Additional hypotheses (SOL), (KER) and (POS): 3TX3l and gOU 

• (nl-KGM) is the nonlinear Klein-Gordon-Maxwell system: ([T]) and m.4\ 

• ^ = (<fi, tp, Aj, Ej) is the dependent variable in the Hamiltonian formulation: 31.41 fand 31.3.11 for zero 
external current case). 
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e electromagnetic coupling constant, 8 external field scaling parameter are both small: ijl.l.3l and £11.2.21 



Scaled external electromagnetic potentials at, and electric and magnetic fields ~E* xt and B^ induced 
5 \ s ) 

Bi JB) 



by external current (p s B , j s B ) : £11.2.21 



^txt — (0) 0, a* 5 , ^txt) represents the external field in the Hamiltonian formulation, and ^ e '*t — 
(0,0,a 5 '* Ef^), its gauge transform by X - SCI and §2~T1 

fwt fu,e are the soliton profile functions in (respectively) the e = case and for non-zero e, while a Uje 
is the Ao component of the electromagnetic potential for soliton solutions: £11.31 

^s,e is the set of Lorentz transformed soliton solutions in Hamiltonian formalism, or ^sc.e in Coulomb 
gauge: §1.1.21 and £11.31 Gauge transformation to Coulomb gauge generated by £: §1.3.41 

A = (A_i, Ao, Ai, . . . , Xq) — (ui, 9, £, u) are parameters for Lorentz transformed solitons: £11.31 

7,^u,Qu,e,9 c ,Z,y (A),iV A ,SandC: EH 

Ostab € K 8 is the stable region of soliton parameter space, where Grillakis-Shatah-Strauss stability 
condition (|39| holds and 5 is positive on the symplectic normal subspace Ny. £11.3.41 

(Ho,0 ), CH,fl) and W^Wn and H^H^ are the symplectic phase spaces and norms: £11.3.1l and £11.2.11 

W, K and 3 are quadratic forms used in stability analysis: £12.31 

W quantity like W but inluding certain nonlinear interaction H parts of the Hamiltonian: £j5l 
7/ oc ,T ,Ti: £EniS2U]and respectively. 



1.2 The External Electromagnetic Field and Scaling 

The external electromagnetic field F^ 1 is induced by the space-time current JJ B = psdt + i%dk according 
to Maxwell's equations, i.e. the first equation of ([T]) with (f> set equal to zero. Introducing an external 
electromagnetic potential, written in lower case symbols, a^dx^ = aodt+ajdx 3 , such that F^ = <9 M a„— c^a^, 
and imposing the Coulomb condition V • a = 0, these equations can be written: 

-Aa =ep B , 

□a =Vd t a - ejs. 

Here pb is the background charge density, j# is the background current density. The associated electric 
field, E ea; t, and magnetic field, B ea;t , are given by 

Eext = ^ - Va , (5) 
at 

&ext = V x a. (6) 
We shall make the following assumptions on the external field: 
(BG) The external electromagnetic potentials are smooth and satisfy: 

max ||V^a M || ioo(K1+3) = L 3 < oo, (7) 

^=0,1,2,3 

(using multi-index notation V" K for arbitrary partial derivatives of order \a\.) 
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It might appear that these assumptions are restrictive: in particular, the assumption that ||a|| Loo ^ R3 xR+ j < 
oo precludes the consideration of a constant magnetic field. However, since we shall scale so that the external 
electric and magnetic fields do not change appreciably over the spread of the soliton, which is exponentially 
localized, these conditions could probably be relaxed with some further work. A more important restriction 
in our study appears to arise in the consideration of the scaling of the the external field, which we discuss 
below, after presenting results on local well-posedness for the (nl-KGM) system in the presence of an external 
field. 



1.2.1 The Cauchy problem for (nl-KGM) in an external field 



Throughout this article we make use of local well-posedness of the (nl-KGM) system in the energy norm. In 
the case that there is no external field and V = this was proved in |13j . In this section we give conditions 
under which this is true in the more general situation of ([1]) considered here. Since our assumptions on the 
external field do not require finite energy it is convenient to subtract off the external field. Thus assume 
given an external electromagnetic potential a$dt + djdxi as above, in Coulomb gauge V • a = 0, which solves 
the inhomogcncous Maxwell equations (TJJ and verifies ([7]) . Write the electromagnetic potential appearing in 
(P) as A M = + Afj,. Then, requiring the Coulomb gauge condition V • A = 0, as is always possible, Q is 
equivalent to the following system: 

4> =ip + ie(a + A )<f), 

tp =A<j) - 2ie(A + a) • - e 2 |A + a| 2 <f> - V'{4>) + ie(a + A )tp, 

□A = (ie<f>, (V - ieA)4>) + VA ~ e 2 |0| 2 a, 
-A-A =(ie4>,ip), 



(8) 



where A = (Ai,A2,A3) is the spatial part of A = A — a. We solve this system in the energy space 
H = H 1 x L 2 x H 1 x L 2 , which is endowed with the energy norm \\^f\\n = IK^^; A, ^)\\h 1 xl 2 xH 1 xL 2 'i see 
^1.1.5l for notation on standard norms. We also define corresponding higher energy norms indexed by s G N 

by 



4,^A,E)|&,= IIV2(^,A,E)|&, 

|a|=0 



(9) 



with corresponding space denoted TC S . We say that the Cauchy problem for (JSJ) is locally well posed in H if 
the following two conditions hold: 

(WP1) given initial data (0(0), ip(0), A(0), A(0)) G H in Coulomb gauge (i.e. divA(O) = 0, divA(O) = 0), 
satisfying 



(#0),V(0),A(0),A(0)) 



H 



< k 



(10) 



there exists T; oc = T; oc (fco) > and a unique solution (((j)(t),tp(t),A(t), A(t)) such that 

(4>{t),tf>(t),A(t),A(t)) G C([0,T loc );H), 
(||DA|| L 2 + Ip^H^)^ < oo. 

(WP2) the solution is continuous with respect to the initial data in that, for another set of initial data 
[(j)l(0),ipl(0), Ai(0), Ai(0)j, which are close in H, and also satisfy (jTUJ) , and the Coulomb gauge 
conditions, the following holds on the common domain of definition [0,T; OC ], for some constant c > 0: 



max 

[o,T ioc ; 



i>x, A- Ai, A- Ai 



< 



H 



(m - mo) MO) - MO), MO) - Ai(o), A(o) - a^o)) 



H 
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As remarked above, in the absence of the external field, and with V = the validity of (WP1)-(WP2) 
was proved in [13]. The general case was addressed in the thesis [17] where it was shown, using in addition 
Strichartz inequalities from [TTJ[33], that (WP1)-(WP2) hold if V is a smooth sub-critical nonlinearity: 

Proposition 1 Suppose V is smooth and that there exists a positive number k S (0, 4) such that, for all 

\V{<j>) - v'Ml <c\<t>-<p\(i + H 4 ~ K + M 4 ~ K ) (li) 

and that V (0) = 0. Assume that the external potential is smooth and verifies ([7]) for every non-negative 
integer j . Then the Cauchy problem for J5]) is well-posed in the sense of (WP1) and (WP2). Further, if the 
initial data lie in TL S for some s > 2 then the solution exists for all time, and remains in TL S , and is smooth 
if the initial data are smooth. 

Remark 2 The Coulomb condition leaves a residual gauge invariance by functions x(A x ) which are har- 
monic in x. (These are either constant or unbounded.) In particular the system ([8]) is invariant under the 
transformation a M i— » a M + 9 M x, (0, tp) <— > e tex (4> 1 ip) if x = a o{t) + ctj(t)x J is linear in x and smooth in 
t. In this case the map ((f), tp) i— » e tx ((/),tp) = is Lipschitz on H 1 x L? . It follows that proposition 

[1] remains valid if the external potential is obtained from one satisfying {7J by gauge transformation by 
X = a (t) + aj(t)x j . 

Remark 3 Notice that when the nonlinearity is determined by a smooth function V whose third derivative 
satisfies: 

£> (3) V(0) < c(l + H 3 ~ K ), for all (12) 

for some c>0, 0<k<3 the conditions of proposition [1] hold, and the Cauchy problem is well-posed. This 
assumption is also sufficient to estimate the nonlinear terms in the perturbation theory developed in fJ2]Sj3] 
and Sj5]of this article. Introduce T((j>) = V'(<p) — m 2 cf> — V[ (4>) = /3(\<f>\)<f> as the nonlinear part of V'(^), with 
V as in the introduction. Then (|12p implies the inequality 

+v)- T'{f)\ < c(l + |/| 3 )(M + \v\ 4 ) where c is a positive constant, (13) 

which is convenient for our use. In fact, for our purposes it would be sufficient to make the following slightly 
more general assumption on T: 

For all / > and for any v, \T'{f + v) - T'{f)\ < c(f- x + f)(\v\ + \v\ 4 ), (14) 

where r, c are positive constants, see [TB]. Of course given a smooth potential V satisfying (Ti"2"]) . let T be as 
just defined, then (TH)) will also hold with r = 1. 

1.2.2 Scaling the external fields 

As already mentioned, we require that the external electric and magnetic fields are approximately constant 
over the soliton. To ensure this, we introduce a scaled version of the external fields. Thus, we have 

4(t, x) = ^a (St, 5x), a s (t, x) = |a(tft, <5x), (15) 

with the scaled external electric and magnetic fields given by: 

Ef xt = E ext (5t, <5x), B 5 ext = B ext (St, <5x). (16) 

Clearly these fields correspond to the following rescaled charge and current densities: 

p s B (t, x) = 5p B (6t, 6x), j s B (t, x) = Sf B (St, Sx). (17) 

Henceforth, we shall almost always refer exclusively to the scaled fields. It remains to choose the length 
scale, 4, over which the external fields change: this is determined by the analysis in Sj5] which bounds the 
deviation of the solution from the modulated soliton. This analysis seems to require two main conditions on 
the scaling of i5 and e: 
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From lemma it seems that we need 



limj=0, (18) 

e^O 



to bound the effect of the scaled external electromagnetic potential. 
• Treatment of the last term in (|137p . seems to suggest that we need 

S 2 

lim — = 0. (19) 

e^O e 

This condition is used to ensure the deviation from the Lorentz force law is small for times of order ^ . 
We will consider the limit e — > with 

5 = e l - k (20) 

for some constant k £ (0, so that both of these conditions hold. It remains to be seen what are the 
optimal conditions for scaling e,S under which the results of this paper hold. 

1.3 Non- Topological Solitons 

We now discuss existence and stability properties of non-topological solitons as solutions of (nl-KGM) in the 
absence of external fields. This means we are here concerned with the (nl-KGM) system with ps = = jb- 
We first discuss the Hamiltonian formulation of (nl-KGM), since that gives the appropriate context in which 
to introduce non-topological solitons. 

1.3.1 Hamiltonian formalism 

It is useful to present the Hamiltonian formalism for the (nl-KGM), not least because it will give us a 
language which we shall use in proving the existence and long-time stability of the non-topological solutions. 
Indeed, as we shall see, from the Hamiltonian point of view, non-topological solitons are relative equilibria, 
and recognizing this fact leads to the identification of the appropriate quantities with which to work. 

In order to define the phase space we recall the standard function spaces defined in §1.1.51 

To start with, consider the nonlinear wave equation in isolation 

D<f> + V'(<f>)=0. (21) 

This can be written as a Hamiltonian system on the phase space Ho = {(</>, 4>) € H 1 x L 2 }, with symplectic 
form f2o((0', tp'), (<t>, ip)) = J (</>' ,ip) — (ip' , (j>) dx, and Hamiltonian 

#oW>,^) = ^/|V0| 2 + 2V(0). (22) 

The corresponding Hamiltonian evolution equations, equivalent to ()21|) . are : 



a I I ) = ( a^,^ )• <:«) 



ip J \ A0-V'(0) 

Next, for (nl-KGM), introduce the phase space 

H = {* = (0, tp, A, E) e H 1 x L 2 x H 1 x L 2 }, (24) 

which is endowed with the norm || T' || = || (</>, ip 1 A, E)|| ff i xL2 x jji xL 2 and the (densely defined, weak) 
symplectic form 

n(* / , *) = J {dp 1 , j>) - W, 0) + A' • E - E' • Adx, (25) 
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where vp' = ((/>', ip', A', E') and similarly for 'J. The (nl-KGM) equations with ps = = js arise formally as 
the Hamiltonian flow on TL associated to the Hamiltonian 



H(4>,1>,A,E) = ± J (|E| 2 + |Vx A| 2 + |^| 2 + |V A ^| 2 + 2V( ( /»)), 
and subject to the constraint: 



C = divE- (ie<f>,ip) = 0. 



(26) 
(27) 



Here Va</> is the covariant derivative of <j> given by Va</> = V</> — ieA(j> and A is the spatial part of the gauge 
field. The equations of motion for the augmented Hamiltonian H\ = H — J AqCq are: 



(28) 



( * \ 




( ip + ieA cj) \ 






A A 0~ V{4>) + ieA iP 


A, 




E, + ViA 


V E * J 




^ A A, - V,(divA) + {iecj>, V A 0) / 



where the "Lagrange multiplier" Aq is identifiable with the temporal part of the gauge field, Aa</> = A</> — 
2ieA • — zediv A(f> + e 2 |A| 2 <f>, i = 1...3, and we have not yet imposed any gauge condition. 



1.3.2 Existence of non-topological solitons: the e = case 

The class of solitary wave solutions of interest is that of non-topological solitons discussed in p~5l Chapter 
7] . These are examples of a special type of solution to a Hamiltonian system with symmetry called relative 
equilibrium: this means that the time evolution is given by an orbit of a one parameter subgroup of the 
symmetry group. For (|23p the Hamiltonian is invariant under the action of S 1 by phase rotation, as long as 
V(4>) = G(\4>\) is a function of \cf>\ only; the charge corresponding to this S 1 action is 

Q(<f>,i>) = J (iil>,4>) dx - 

A relative equilibrium is then a solution of the form (<fi, ip) = e lu)t (f u] {x),iujf UJ ) 1 where f u is a real valued 
function which satisfies an elliptic equation. These solutions are critical points of the functional Hq + ujQ, 
often called the augmented Hamiltonian in this context. We consider G of the form 

777 2 

G(.f) = —f - U(f), With U(f) = tp(t)dt. 

then the equation satisfied by f u is — A/ w + (m 2 — u) 2 )f u = /3(/^)/w This equation typically has many 
solutions (see [3] and references therein), but we are only interested in positive, radially symmetric solutions 
because it is these which are dynamically stable: these are sometimes called the ground state solitons. Thus, 
crucial to our analysis is the following hypothesis on existence and uniqueness of the e = ground state 
soliton: 

(SOL) For u! 2 < to 2 , there exists a unique positive radial function f u e H 4 (R 3 ) which solves 

(-A + to 2 -^ 2 )/ w = /3(/ w )/ w . 

Theorem 4 The existence part of (SOL) holds under the following conditions: 

U'(f) = -U'(-f) and U' E C 1 (M)nC 2 ((0,oo)), (29) 
U'(0) = U"(0) = and 3s € (0, 1) : lim f s U"'(f) = 0, (30) 
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3C>0:t/(C)>^ T ^C 2 , (31) 
U'(f) 

lim = 0. (32 

The uniqueness part of (SOL) holds under the additional conditions: 

(Ul) 3h>0:0< f <h=^ U'(f)< (to 2 - cu 2 )f 

and h < f < oo => U'(f) > (to 2 - u?)f 
and U"(h) - (to 2 - w 2 ) > 0, 

and that 

(U2) Forh > h,3X = \(h) G C[(Zi, oo), R+] 

suc/i i/iat 2(m 2 - w 2 )/ + XfU'(f) — (A + 2)J7'(/) 

«s non-negative on (0,^) and non-positive on (£2,00). 

Proof The existence part of this hypothesis was proved in [3] under the given conditions on the nonlinearity. 
It was shown in several articles (see for example, [18], where further references are given), that these solutions 
are unique under the given additional conditions. □ 
The following two operators, L±(uS), which appear on linearizing (|23p about the soliton solution, are 
crucial to an understanding of the stability and dynamical properties of the e = soliton: 



(33) 



L+H = -A + m 2 -cj 2 - p(f u ) - p(J u )f u , 
L-(u) = -A + to 2 - lj 2 - (3(f u ). 

We make the following hypothesis on L + (oj): 

(KER) The kernel of L+(ui) is empty in H 2 (R 3 ). 

(Recall that H* was defined as the space of radial Sobolev H s functions, immediately following Q.) 
Theorem 5 The hypothesis (KER) is valid under the conditions (U1)-(U2). 

Proof See |18) : establishing (KER) is a crucial step in proving uniqueness of the positive function f u . □ 
The operators L± also determine stability properties of the soliton. For proving stability the following 
spectral assumption is used: 

(SI) The subspace in which L + is strictly negative is one dimensional, 

This assumption is valid for the ground state solitons f u obtained by the constrained minimisation 
technique of [3], because they are minimizers subject to a single constraint, see |24j (where a direct proof in 
the pure power case is also given). 

Some additional more technical results on the solitons can be found in appendix lA.il 



1.3.3 Existence of non-topological solitons: the general case 

We now show that for small values of the coupling constant e the ground state solitons just discussed can be 
continued (via the implicit function theorem) to give soliton solutions of (|28p . The properties of the e = 
soliton needed to achieve this were detailed already in i ll. 3. 21 As shown in [TJ [2] it is also possible to obtain 
soliton solutions for systems like (|28|) by variational techniques applied within the class of radial functions, 
but for present purposes we prefer to use the implicit function theorem so that we can carry over stability 
information from the e = case, which seems to be hard to obtain otherwise. 
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Generalizing the class of non-topological solitons to the case of the gauge invariant system (|2"5)) leads us 
to search for solutions to ([28]) of the form 





* \ 




( 


Exp[uot]f u>e \ 










Exp[iu;t]i(uj - ea Ut ^)f u>e 




A 









V 


E J 




\ 





(34) 



where we have emphasized the dependence on the parameters to and e; we will assume the functions f u<e 
and otu, e to be radially symmetric. It can easily be checked that this gives a solution to (I28p with Aq = a^.e 
as long as the functions f u>e and a u . e satisfy 

-Aa u , e + e 2 /^ e a w , e - etof% <e = 0, (35) 
-Af u , e - U'{f u , e ) + (m 2 - (w + ca,, e ) 2 )/,, e - 0. (36) 

The first of these equations implies Co = 0. It can readily be checked that if a gauge transformation is 
made to put the solution thus obtained into temporal gauge, Aq = 0, then its time dependence amounts to 
the action of the one parameter group of gauge transformations e l ( u ~ eoLu, - e * >t , so that it is indeed a relative 
equilibrium solution as defined above. 

Theorem 6 Q17J) Assume that the hypotheses (SOL) and (KER) hold for to q with ujq 2 < m 2 . Then, there 
exists a neighbourhood U of ujq such that for to G U , there is a number e(to) > such that for to G U, |e| < e(to), 
there exists / w e G H 2 (M. 3 ) such that 

- Mu,e + m 2 L,e - (W - ea^e) 2 f Ut e = /3(f u ,e)fu,e, (37) 

where a U)e G H^:(M. 3 ) is a non-local function of f u>e uniquely determined by 

-Aa^ + e^A^we/^- (38) 
In addition the map lu i— * f ue is C 2 from U to H 2 . 

1.3.4 Stability in the absence of an external field 

The stability of the solutions to (|2"3"]) of the form e luJt f LU (x) was first considered in [5D1[TD] where it was proved 
that the positive radial solution was stable, with respect to radially symmetric perturbations of the initial 
data, as long as 

a w (^||/ w || 2 2 ) <0. (39) 

It was also shown that the solutions are unstable when this quantity is positive. In |24j an alternative, modu- 
lational, approach to stability was adopted along the lines of [27] , with the aim, both of generalizing previous 
stability results to prove stability of uniformly moving solutions with respect to arbitrary (non-symmetric) 
perturbations, and also of providing techniques which could provide useful information in dynamically non- 
trivial settings. The presence of external fields is an example of the latter circumstance, and so the analysis 
in this article is based on that in [24] . which we will now summarize. It turns out that the condition ()39() 
implies the strict positivity of the Hessian of the augmented Hamiltonian on the symplectic normal space to 
the space of solitons. To explain this properly in the generality needed it is necessary to consider the action 
of the Poincare (or inhomogeneous Lorentz) group 

Action of the Poincare group on the solitons. The equations (|28p are Poincare covariant. The action of the 
Poincare group on the radial soliton (|34[) gives a family of functions depending smoothly on eight parameters 
{AA}i = _ l5 with 

A = (A_i,A 0l Ai,... ) A 6 ) = (w,0,£,u) (40) 
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determining (respectively) the frequency, the phase, the centre and the velocity of the soliton. Explicitly: 



*s, e (x; A) 



/ Ex P [ie](f u , e (z)) \ 

Ex P [ie}(ij{LJ - ea w , B (Z))/ w , e (Z) - 7 u • V z /,, e (Z)) 
-7ua Wie (Z) 
-(ifu + 7Qu)V z aa,, e (Z) 



Here the projection operators P u : M 3 
and 



and Q u 



are defined by (P u )i 



J 



\u\- 



(41) 



and Q u = 1-P U , 



Z(x,A)=7P u (x-€) + Q u (x-€), 



(42) 



6(x,A) = 6»-cju-Z, 

with 7(11) = -, 1 The parameters are required to lie in the set Oct 8 defined by 

6 = {(u,9,€,u) C R 8 : |u| < 1 and uj 2 < m 2 }. 
The parameter range corresponding to stable solitons is 

Ostab = {{uj, 6, £, u) C O : condition (J39J) holds}. 



(43) 



(44) 



(45) 



The Poincare covariance of the equations of motion ([28)) implies that the solitons given by (|4Tj) form an eight 
parameter family of solutions t 1— ► 4 , 5. e (x; A(i)) of (j2"5| as long as ^A = Vo(A), where Vq is the vector field 

on O defined by 

K,(A) = (0,-,u,0), (46) 

for A = (ui, 8, £, u). 

The case of the nonlinear wave equation (|23[) can be obtained by putting e = in the first two components 
of the formulae just given. Simplifying to this case we obtain an eight parameter family of functions, 



,o,^ s ,o)(x; A) = e ie (UZ), (i>yuf u (Z) - 7 u • V Z / W (Z))) 



(47) 



such that 

solves ([23)) . as long as ^A = Vq(A), with Vq as above. 

Stability for e — (nonlinear Klein- Gordon). The starting point for stability analysis is the observation that 
(</>s,0j V'S.o) is a critical point of the augmented Hamiltonian 



F (cj>, i>- A) = H (<f>, i>) + mTT(0, 4,) + -Q(4>, 

1 



(48) 



where Hq,Q are the functionals defined above, and IT are the momenta IT(</>, ip) = f(ip,di(j)}dx. The 
Hessian of Fq at (0s, 0, V'S.o) is a quadratic form depending upon A: 

3(0, fa A) = D 2 F^ Sfil ^s,o; A)((0, (0, i>)). 

Introduce the subspacc 

N x = e H 1 x L 2 : n ((^^),9A(^,o,^,o)(A)) = 0} (49) 

then the following hypothesis is crucial for stability: 
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(POS) For each compact K C O s tab 3t* = (K) > such that 
E(4>,i>;X)>T4(^i>)\\ 2 HlxL2 for all e N x . 



Remark 7 O s t a b is the set of parameter values corresponding to stable solitons, which are obtained as 
Poincare transforms of solitons e tu}t : / w with oj such that ([39| holds. 



Theorem 8 (|24|) If the nonlinearity satisfies the conditions given in Q1.3.2\ then (POS) is true. Further- 
more, solitons of (I23p corresponding to frequencies oj such that (|39j) holds are modulationally stable with 
respect to small, arbitrary perturbation in energy norm. To be precise, consider the initial value problem for 
(f2"3")l with initial data close to a soliton (<fis,o, V>S,o)( - ; A(0)) with A(0) £ O s t a b, in the sense that 

e = || (0(0, •), ^(0, 0) " (0s,o, ^s,o)(-; A(0)) || Wq 
is sufficiently small. Then there exists a global solution which satisfies: 

sup 1 1 (0(t, •), VK*, ')) - (0s,o, ^s,o)(s A(i)) || Wq < ce, 



(50) 



/or some C 1 curve t i— > A(i) € O stab- 
Stability for small e (nonlinear Klein- Gordon- Maxwell). It was shown in [17j . that stability holds also for 
solitons in ([28]) under the condition (|39|) , for sufficiently small values of the electromagnetic coupling constant 
e. This was proved using the Coulomb condition, so we first write down the soliton solutions (|41[) in Coulomb 
gauge. (The Coulomb condition is not invariant under Lorentz boosts, therefore, it is necessary to perform a 
gauge transformation to move the Lorentz boosted solitons into the Coulomb gauge). The Lorentz boosted 
solitons \& sc,e in the Coulomb gauge have the form 



/ 


0SC,e(x) ^ 


( 










AsCe(x) 




\ 


EsC,e(x) J 


V 



Exp{iQcKUA z )) 
Exp[iQ c ]{il{u - ea w , e (Z))/ w , e (Z) - 7 u.V z / u , e (Z)) 
-7ua Wie (Z) + VC 
-(^u+7Qu)V z a Wje (Z) 



\ 



(51) 



where ®c = B + ie(, and £(x; A) is a solution of 

- AC= -7u-Va w , e (Z). 



(52) 



It is a smooth function of x and also depends smoothly on A; requiring that V£ £ i p ,p > 3 fixes it up to a 
constant. Some estimates for £ are given in appendix IA. 1.21 The temporal part of the gauge field is given 

by 

(Asc, e )o = ja u , e (Z) + C - ja u , e (Z) + Vb(A) • d x (. 

Theorem 9 Q17J) In the situation of the previous theorem the solitons (1511) of (|28[) corresponding to fre- 
quencies oj such that (|39[) holds are, for sufficiently small \e\, modulationally stable in Coulomb gauge with 
respect to small, arbitrary perturbation of the initial data in energy norm \\ ■ defined in (|24p . The stability 
is in the same sense as in the previous theorem, see \17^ for full details. 
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1.4 The main theorems 



We now explain and state our main results on the interaction of the solitons of £11.31 with the scaled external 
electromagnetic field of HI. 21 We write the total electromagnetic potential as A = A^dx 11 (as described 
in i jl.l.ip with corresponding electric field Ej = d t Aj — djA . The potential A will be formed from three 
constituents: 

1. the external field, produced by a background charge p B and current j^, and scaled as described in 
CCS 



2. the soliton contribution, as described in £|1.3l but with parameters X(t) varying in a dynamically deter- 
mined way, 

3. an additional component produced by interaction of the initial data with the two previous components. 
This component is not explicitly given, and must be estimated. 

Similarly, the solitonic field will be made up of a component which is the moving soliton, and a remainder 
produced by interactions, which must be estimated. 

It is convenient to write the (nl-KGM) equations in first order form. Including the scaled background 
current density, the equations read: 



/ <P \ 



di 



tp + ieAo</> 
A A 0-V'(4>)+ieAoV> 
E» + ViAo 
\ AA 4 + (ie4>, W A (f>) - ejfj / 



A t 
V ) 

with the Coulomb gauge condition imposed. These equations are to be solved with the Gauss law 

divE — (ie<f>, ip) = p s B , 

as a constraint. We shall abbreviate a general solution by making use of the following definition: 

* = (0, ^, Aj,Ej) , 

with i £ {1,2,3}. 

Using this Hamiltonian formulation with \l/ as dynamical variable we write the external field 



(53) 



(54) 



(55) 



(o,oy,EL t ). 



It will be convenient also to have the freedom of applying a gauge transformation x(i, x) to this: 



(0,0,a^,E^ t ), 



with a S:X = a s + A\- The aim is now to construct a solution ^ to (|53|) consisting of with a soliton 
^SC,e(A) superimposed. We choose the gauge transformation \ so that the transformed external electro- 
magnetic potentials vanish along the world-line of the soliton x = £(t); in particular, at t — we will 
choose 

X (0, x) = xo(x) = -(x - {(0)) • a 5 (0, 
so that a 5 >x°(0,x) = a 5 (0, x) - a 5 (0, C(0)). 



1.4.1 Stability in the presence of an external field 

The following theorem asserts the long time stability, under the influence of an external field, of stable 
solitons to (|53|) . Recall that the stable solitons are those parametrised by A € O s tab, so that |39|) and hence 
(POS) hold, and they are stable by theorems [5] and [5] in the absence of an external field. 
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Theorem 10 Assume that the nonlinearity satisfies the hypotheses (Hl)-(HS), and also is such that the 
hypotheses (SOL), (KER) and (POS) in ^1.3.2\ hold. In addition, assume that the external field satisfies the 
assumptions in JHH Suppose further that the scaling parameters satisfy S 2 = o(e), e = o(l) and e = 0(6). 

(i) Consider initial data of the form 

*(0) = (^ sc ,e(A(0)),^sc,e(A(0)),A i (0),E i (0)) 

where A(0) G O s t a b corresponds to a stable soliton (which verifies (POS)). It follows that, if e is sufficiently 
small and 

||*(0) - *£?(0) - * sc ,e(A(0))||* = o(e), (56) 

there exists 

• a positive number Tq > 0, independent of e, 

• a C 1 gauge transformation \(t,x) defined in (|63|) . linear in x af eac/i iime t, satisfying x(0,x) = xo( x ) 

• a curue A(i) € ^([O, S], O sto6 ), and 

• a distributional solution ^(t) of H53\) . 
such that 

*(t)-*^(t)eC([0, 

|e| 

and 

sup ||*(t)-^(t)-*so, e (A(t))||^ = o(e), (57) 
te[o,^f] 

Furthermore, X(t) satisfies a system of ordinary differential equations given by 1116\) with \dtX — Vq(A)| = 
0{e). The time component of the potential Aq is determined by the Coulomb condition and the Gauss law, 
(jTTj) , and has properties detailed in SJH 

(ii) More generally, the same conclusions hold for initial data sufficiently close to a stable soliton in an 
appropriate sense: see \ \2.4-3\ for a precise statement. 

This theorem is proved in Sj2] 

Remark 11 As explained in 11.3.21 if the nonlinear potential satisfies (|2"9]) -([32" j) . U(l), U(2), 5(1) above in 
addition to (H1)-(H3) then the conditions (SOL), (KER) and (POS) all hold. 

1.4.2 Motion in the presence of an external field: the Lorentz force 

The previous theorem provides ordinary differential equations (| 1 which determine the evolution of the 
soliton parameters. A detailed investigation of these equations allows us to deduce an equation of motion 
for the soliton, which is expected to be the Lorentz force law for a moving charge, at least to highest order 
in e. As remarked earlier, if the analysis were carried out explicitly to higher order in e, corrections would 
be expected to appear, in particular due to the back reaction of the soliton's electromagnetic field on itself. 
However, these are not expected to appear in the 0{e) force law, and the following theorem validates this: 

Theorem 12 Assume the hypotheses and conclusions of theorem rjTJ hold, and let A = (u>,6, £, u) be the 
parameters of the soliton v I'sc,e(A). Then, on the interval [0, Sj-], the centre and velocity of the soliton evolve 
according to the equations: 

u + o(e) (58) 
eQ s (EL t (*,0 + u x B s ext (t,$) + o(e), (59) 



- (M s7 (u)u) 
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where the mass of the soliton, Mg, is given by 

Ms = l\\Vfu,t L2 + 0J 2 \\fJ 2 L2> (60) 
and the charge of the soliton is given by 

Qs = J (w-ea)/* e . (61) 

This theorem is proved in <|4] 

Remark 13 Observe that, since we have scaled the external field so that E s ext and Bf It are independent of 
e, the soliton undergoes O(l) motion on the time interval [0, jk] according to the Lorentz force law. 



2 Stability: proof of theorem [10 

In this section we explain the proof of theorem 1 101 making use of results which are proved separately in 
$3] and Sj5] Throughout this section the hypotheses of theorem [10] are understood to hold without explicit 
mention. Also we may assume, without loss of generality, that the solution is smooth in the course of the 
following calculations: since finite energy solutions can be approximated by smooth ones by (WP2) in £H.2.1| 
and all the bounds we use depend only on the energy norm, this implies the result for finite energy initial 
data as in theorem [T0l 



2.1 Beginning of proof of theorem 1101 
2.1.1 Ansatz for the solution 

We make an ansatz for a solution = ^ S e ' xt (t) + ^> sc,e(A(t)) + Perturbation, which is close to a soliton 
with time varying (modulating) parameters A(t), in the background external field \l/f^(t). Explicitly the 
ansatz reads: 



0sc,e (A(t)) + Exp[iQ c ]v 
^sc.e (A(£)) + Exp[iQ c ]w 
(A sae )^\(t))+a^ + A 
V (E SC ,e) 3 (A(t)) + (E s ext ) 3 

Notice that we have included here an ansatz for the temporal part of the potential 



/ 0(*,x) \ / 
^(i,x) 
A M (<,x) 



1 m 
E 3 ) 



(62) 



Since we have imposed 



the Coulomb gauge throughout, it follows that div A = 0. The choice of the gauge transformation \ is: 

X (t, x) = -(x - i) ■ a 5 (t, £) - ( a s (s, {(a)) + £(s) ■ a 5 ( S , £)ds. (63) 



This is chosen so that the gauge transformed external potentials vanish along the world line of the soliton: 



These imply 



4< x (t, x) = a 6 (t, x) - a a Q (t, {(t)) -(x- • a" (t, €(t)), 



(64) 



(65) 



exhibiting the claimed vanishing of a*' x along the soliton's world line. This allows certain quantities to be 
proved to be bounded in the course of the proof. Notice that \ IS linear (and so harmonic) in x, and so 
preserves the Coulomb condition (see remark [2]). 
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There is clearly a redundancy in our ansatz, in that X(t) is so far completely undetermined. The appro- 
priate choice of X(t) is dictated by the requirement that the solution be close to a soliton determined by the 
parameters A(t), i.e. by the requirement that we have good bounds for field perturbation (v, w, A, E). This 
is carried out in with the main results summarized next in Q2.2\ First we write explicitly the equations 
for the (v, w, A, E), and give some bounds for the inhomogeneous terms in these equations. 

2.1.2 Equations for the perturbations of the fields 

dtv + i (uj + h) v = w + ji, (66) 

d t w + i(oj~/ + h)w = —M\v + (67) 

d t A = E + j 3 , (68) 

<9 t E = AA + j 4 , (69) 
where the inhomogeneous terms h,j\,...,j4 and J\f are defined in ij2.1.31 and M\ is the operator 

M x v = (-A a + m 2 + 7 V>|> + 2zc^ 7 u • V x v - 0(f u )v - Uf3'(U)Viv. (70) 

The last two terms have been chosen to depend on the e = profile function f u , rather than f ue , so that 
it is possible to make direct use of the stability assumption (POS) in HI. 3. 41 (This choice is reflected in the 
expression for the inhomogeneous term J\f in (|76p and its corresponding estimate in (|95p). 

In addition to these evolution equations, the fields are constrained to satisfy the Gauss law (f27|) . which 
takes the form: 

divE = -AA = e (iExp[-i®c](f>SC,e-,w) + e (iv, Exp[-~iQ c ]^SC,e + w) . (71) 

Under finite energy assumptions this equation has a unique solution with Aq £ H , this defines uniquely A$ 
as a nonlocal function of v, w, X at each time. Estimates for Ao are given in lemma [39] 

2.1.3 Inhomogeneous terms in the field perturbation equations (|66]) - (j69|) 

The following quantity appears in both (|6"6")) and (p71) : 

h = 6 C - W7 - e(A S c,e)o - ea S Q ' x - eA Q . (72) 
The inhomogeneous term in (|66p is ji = j[ + j[ ! + , where 

A = -(A-V o (A))-e- ie ^ A 0sc,e (73) 

A 1 = (74) 

A = ieAof u , e . (75) 
The inhomogeneous terms in (|67p are 

M(f u ,e,U,v) = 0{\f«,* + v\){U,e + v)-(3{\U, e \)U,e (76) 

-P(\f u \)v-f u 0'(\f u \nv, 

and j 2 = jl + 3" + J2" + A 2 V + J°2 where 

A = -^-V (X))-e- l ^d^sc.e 

AA = eR/ W:e + iea^ e - ie ^sc,e, (77) 
3A 1 = eRu + Sv (78) 
A 2 V = e- ie ° [A A (e ie ^) + (A A - A Asc> J0 sc , e ] - AA ~ A" (79) 
3° 2 = ieA e~ ie ^sc,e- (80) 
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Here, the operators R, S are given by 



Rv = 2i(a S ' x ) .(i-y(u-ea u , e )u- V)v - e\a. S ' x \ 2 v, (81) 
St; = 2iea u . e '-fU ■ Vf + iej (u • Va^.e) v + 2ej 2 \u\ 2 aja UJ . e v — e 2 (7a Wie |u|) 2 v. 
(In verifying these formulas, it is helpful to note that by the exact solutions in £11.3.41 

e" l9c (V - ieA sc , e - iea StX )e i@e v = Vu - i(-y(w - ea w , e )u + ea S ' x )v, 

and a similar formula for the second derivatives.) 

The inhomogeneous term in (|68[) is J3 = J3 + J3 where 

3 I 3 = -(\-V (\))-d X A SC ,e, 

Jl = Vl - 

and in ([69]) we have j 4 = j{ + ]" + j{ n + j{ V + 3°, with ]\ = and 

.A = -(\-V (X))-dxE sc , e , (83) 

jl 1 = -e 2 |^, e |V>* (84) 

j{" = e^SVAfc + ^^-e'l/./A, (85) 

jf = e(^ sc , e ,V A (e 40 ^)). (86) 

To clarify the structure of these terms it is helpful to insert the ansatz (|6"2"j) into the Hamiltonian (f2"6")) 
and write H — f V = X) ra =o , where f/^™) has homogeneity n in (u, A). (The terms of degree larger than 
two arise solely from i J |Va</>| 2 -) Then the pieces of (resp. J4), which are of degree n € {1, 2, 3} in (v, A) 
arise, respectively, as the Frechet derivatives —D v H^ n+1 \ (resp. — D A H^ n+1 ^). The nonlinear potential V 
only appears through M\v and A/" in (|67l) . With this understood we now introduce notation for the various 
terms arising in (f6"T)l and (|6"9"]) , organized according to their homogeneity. Let H = Ylti=2 > then in (f6"T|) 
the corresponding terms are 

-D v H = -M x v+jI ii +j i 2 v , 

where 

M^w = (-A a + j 2 uj 2 \u\ 2 )v + 2iw 7 u • V x v = -e- m A(e m v) (87) 

with Q the soliton phase factor in (j4"3"|) . Notice that the operator M° consists of those terms in (|70l) . which 
do not arise from the V term in the energy, because we have so far excluded this term in our expansion 
(which is of H — J V). However, it is convenient to put back in the quadratic parts of the Taylor expansion 
of V, but expanded around f u (the e = soliton), so as to obtain the M\ operator which appears in l|67[) . 
Thus we let 

H = H + l -D 2 V{f^){v, v)=H+ 1 - [m 2 H 2 - 0(J u )\v\ a - f u p(f u )($v) 

so that, using the same notation for the homogeneous components of H as for H, we have H^ n ' = for 
n > 2 and — = ^D v [f u )[v, v). In (f6"9")l the corresponding terms are 

-D A H = AA + jl n +J 1 /. 
To write these terms explicitly we introduce a multilinear notation as follows. 

{-D v H,-D A H)=M^{v,A)+B^(v,A)+M^(v,A), (88) 

where M^(v,A) is a homogeneous degree n function of (v,A), as indicated by the superscript. We will 
define 1^ : fl^R^C) 8 ij\(M 3 ;R 3 ) >-> H^iR 3 ^) ij-^R 3 ; R 3 ), where by ff" 1 (resp. ff" 1 ) we mean 
the dual space of H 1 (resp. H .) Explicitly: 

B (1) (v, A) = (B u t. + B12 A, M 21 v + B 22 A), (89) 



(82) 
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where 

Muv = —M\v + eKv + Sv, (90) 
and the operators R and S are as just defined. Next 

B 12 A = -2ef u , e (7 (w - ea w , e ) u + ea S ' x ) • A - 2zeA • V/ W , e , 

B 2 i v = -2e/ w , e (7 (w - ea w , e ) u + ea" 5 ^) + (ieu, V/ W , e ) + (ie/ w>e , Vu), 

and finally, B 22 A = AA — e 2 / w-e A. Since divA = integration by parts yields (A,B 2 i v) L 2 = (u,Bi2 A) L 2, 
and 



£(2) 



--((,,A),B (1) KA); iL 



(91) 



iVAf + e 2 |/ w , e A| 2 + |V« - t(7 (w - ea w , e ) u + ea a ' x ) W 

+ i(m 2 |«| 2 - 0{f u )\v\ 2 - WiUWvf) 

- 2{ieAf< JJ ,e, V« — i(7 (w — eow) u + ea 5 ' x )u) 

- 2(ieAu, V/ W , e - i(7 (w - ea w , e ) u + ea" 5 '*)/^) 



Next, the quadratic terms in the equations can be expressed in terms of a rank three symmetric tensor 
B (2) : (H 1 (R 3 ;C) ® H 1 (M 3 ;R 3 )) ff- 1 (R 3 ;C) © H' 1 (K 3 ;M 3 ) 
which is given explicitly by 

B ( 2 ) f fa = ( B iii[ w ^]+ B ii2b ; A ]+ B i2i[A,«]+B 122 [A,A], 
lU ' j V B2n[«,«] +B 2 i 2 b,A] +B 22 i[A,u] +B 222 [A,A] 



where B m 



^222 = 0, and 

B112 [v, A] = —ev (7 (no — ea^.e) u + e(a <5,x )) • A - ieVv ■ A, 

B121 [A, v] = —ev (7 (lo — e« W) e) u + e(a <5 ' x )) • A - ieVv ■ A, 
j [A, A] = -e 2 /u,,e|A| 2 , 



»122[ 



and 

along with 



8211 [u, u] = -e (7 (to - ea U!e ) u + e(a ,x )) H 2 + (iew, Vu) 



B 221 [A,i>] =B 212 [«,A] - -e 2 (/„,»,«) A. 
These terms are obtained by differentiation of the cubic part of the expanded Hamiltonian, which is 

#(3) = 



--((«,A),BW(t),A); y 
= (Vw — 1711(0; — ea u , e )v — iea. Sa v, —ieAv) L 2 + e 2 (Af^, e , Av) 1?. 
Finally the cubic terms in the equations arise by differentiation of the quartic part of the Hamiltonian 



# = 1 ((„,a),i 13 '(«,a); [! 



1 



y / |A| 2 M 2 , 
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and are determined by a rank four tensor, 

B (3) : (if 1 (R 3 ; C) H l (R 3 ; R 3 )^ t-y if" 1 (R 3 ; C) if" 1 (R 3 ; R 3 ) 

which, using an identical notation to the rank three case, has as its only non-zero entries 

-e 2 

B1122 = — (92) 
and the other entries obtained by permuting the indices. 

2.1.4 Some bounds for the inhomogeneous terms 

We record here some simple bounds for the quantities defined above: 

• WWlp + WtfWis = 0(e) and ||jf \\ LP = 0(e 2 ) for every p £ [1, oo] by (HHZD , dHD , 

• ||fr/w,e||Lp = 0(e + |A — Vb| + e|Ao|L9), for any q > 3, which can be read off from (|72|) . using results 
from appendices I A. 1.21 A. 2. 21 and IA. 2. 11 and the assumptions on the applied fields. A can be bounded 
in L q , q > 3 by appendix [A. 2. 2 1 

• It is possible to write h = hi — eAg with || V/ii||l°° = 0(e+| A— Vb|) and V A$ bounded in L p , p £ (3/2,3], 
by appendix IA. 2. 21 

Finally, consider TV: by lemma [35] we can write 

M(U, B ,f U >v) = P(\f u , e + V\)(f u , e + V)-I3(\f u , e \)f u , e 

-P(\f u ,e\)v ~ U,e/3'(\fu>,e\nV + 0(e 2 \v\) 

= M{L,eiL^v)+0{e 2 \v\) (93) 
Using the condition (fT2")l . or more generally (fT"3"|) , and the fundamental theorem of calculus, we can estimate 

Wif, f,v)\<c(l + |/| 3 ) (\v\ 2 + \v\ 5 ) , (94) 
for any /. Therefore, choosing / = / w>e , which is bounded, and using (|93[) we have 

\mUe,fu,v)\< Cl (\v\ 2 + \v\ 5 )+c 2 e 2 \v\. (95) 

2.2 Results from modulation theory 

The assumptions on the nonlinearity under which we are working ensure that the Cauchy problem for (|53[) 
is locally well-posed in the sense of (WP1) and (WP2), see £11.2.11 Since so far \ is unknown (since X(t) 
and hence £(t) are not yet determined) we cannot solve directly for 'J = (</>, ip, Aj, Kj) in the background 
potential a s ^ x . Instead we exploit gauge invariance and solve for 

* = (fafak^Ej) = (e-' iex </),e- iex V, Aj - djx,®i) = e~ iex ■ (96) 

in the potential a*, which is known. (Since x(f,x) is harmonic in x this gauge transformation preserves 
both the equations ([53]) and the Coulomb gauge condition (see remark [2])). By proposition [T] on local 
well-posedness, there exists a time Ti oc > and unique solution to (|28p ) with 

(*-** It )eC([0,r te ];?i), (97) 

with initial data 

*(0) = (e- <e *V(fJ),e- <e *°V>(0), A .( ) _ fi^E^O)). (98) 
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Once \(t) = (u)(t),6(t),£(t),u(t)), and hence x(*)i is determined, then *S?(t) is obtained from ^(t) by the 
above relation. As remarked previously, by proposition [T] these solutions can be approximated in energy 
norm by smooth solutions evolving in any of the spaces Ti. s of ^ (after subtracting off the background 
field). Thus, although the statement and proof of theorem 1101 involve only the energy norm, it is permissible 
to assume smoothness of the solutions throughout the proof. 

We now state a theorem which asserts that it is possible to choose the soliton parameters X(t) in such 
a way that the quantity W defined in p03[) is equivalent to the energy norm. This is achieved by choosing 
X(t) in such a way that the pair (v, w) satisfies some conditions which are equivalent to those in (|4*^1) (after 
adjusting the phase). 

Theorem 14 (a) Let ^ be a solution to the Cauchy problem for \5S\) satisfying (I97p with initial data 
(|98p with \l/(0) as described in theorem \1(A Then, for sufficiently small e, there exists T\ > and A G 
C 1 ([0, Ti]; Ostab) with the following properties. On the interval [0, T{\ define ^(t) — ((f)(t),ip(t),Aj(t),W,j(t)) 
by (1631) and (|96[) . Then it is possible to write $ in the form (|62[) where v,w are constrained to satisfy 

n ((v,w),^s,o^s,o)) =0, (99) 

where we define 

dx(jy Sfl = Exp[~i (Q)]dx4>s,o, (100) 

and likewise for d\tps,o- Furthermore, the function t i— » A(t) solves a system of differential equations (|116p . 
The condition (|99|) is equivalent to requiring (</> — 4>sc,e,ip — %j)sc,e) £ N\. 
(b) If e and \\ (v,w, A, E)||-^ are sufficiently small, then 

|A-V (A)| = 0(e + \\(v, w, A, E)f H ) , (101) 

so that, in particular, if \\(v, w, A, E)||^ = 0(e) then 

\X-V (A)\ = O(e). (102) 
Proof This is a consequence of the lemmas in $3l □ 

2.3 The main growth estimate 

As discussed in §1.3. 4[ the natural quantity for stability and perturbation analyses of the solitons (l5T|) is the 
Hessian of the augmented Hamiltonian. Here we modify this quantity to take account of the phase shifts in 
(|62|) . and discard terms which are formally O(e), leading us to the introduction of the following quadratic 
form: 

W(v, w, A, E; A) = K + E, (103) 

where 

K(A, E; A) = \ (||E||i 2 + ||V x A||| 2 + 2(E, (u • V) A) L .) , (104) 

and 

E(v, w; A) = - (\\w - i!UJv\\ 2 L2 + (v, M\ - i 2 oj 2 )v) l2 + 2 (w, u • V«) L2 ) , (105) 
where M\ is as defined in (|7D| . 

Theorem 15 (Equivalence of W and energy norm) Suppose that the nonlinearity is such that (Hl)- 
(H3) and (SOL), (KER) and (POS) hold. Suppose further that A lies in a compact subset, K, of O s tab- 
Then the quadratic form W just defined, is equivalent uniformly on K to \\(v, w, A, E)||^ provided that 
(v,w) satisfy the constraints (|99p . 
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Proof This is essentially theorem 2.7 in [2"1] . Since there is no coupling in W between (v, w) and (A, E), it is 
only necessary to show separately the equivalence of S and K to the corresponding parts of ||(i>, w, A, E)||^. 
For K this can be achieved by completing the square (since || V x A|| i 2 = || A||^-i by the Coulomb condition), 
while for S it is an immediate consequence of (POS). □ 

Theorem 16 (Main growth estimate) Assume given a solution to the Cauchy problem for (28\) for which 
theorem \lJ\ applies on an interval [0, S-] for some fixed positive T2. Assume that \(t) S K, a compact subset 

of O s tab, so that by theorem \15\ there exists c\ > such that, 

-W<\\(v,w,A,E)\\n<ciW, (106) 

Cl 

on [0, S]. Assume further that there exist C2 > 0, C3 > such that that S 2 < C2 |e| and and W < C3 1 e 1 7 and 
that e = o(5). It follows that, for sufficiently small e, there exists C4 > such that, on [0, t4] 

W(t) < c 4 (VF(0) + e 2 + 5 2 ) exp(c 4 |e|i). (107) 
Proof See S0 □ 



2.4 Completion of the proof of theorem [TUl 
2.4.1 Local solution verifying constraints 

For simplicity of exposition we first prove part (i) of the theorem, i.e. we consider initial data ^(O) consisting 
of an exact soliton as in (f5Tj) determined by parameters A(0) = (0(0), w(0), u(0), £(0)) € O s tab, with oj(0) 
satisfying the stability condition. On account of the applied fields there will be a non-trivial evolution starting 
from this initial value. Applying the local existence theorem [TJ and theorem [TJ] as in W2.21 we deduce the 
existence a positive time T\ > such that on the interval [0, Ti] there is a solution to the Cauchy problem 
which can be written as in ([62]) where v(0) = = w(0), and (v(t),w(t)) satisfy the constraints (|109[) (or 
(fM]0 , and t 1— > A(t) solves (|116[) . We may assume that X(t) € K, a fixed compact subset of O s tab, so that 
(|T06l) holds. 



2.4.2 Growth of the energy norm 

Since we have a local solution satisfying the constraints we can assume that the conclusions of theorem 
1151 hold. Furthermore, by continuity we may assume (making T\ smaller if need be) that on this interval 
W(t) < C3 1 e I , and (|106[) holds. Now apply the growth estimate in theorem [TBI 

W(t) < c 4 (W(0) + e 2 + S 2 ) exp(c 4 |e|i), 

to deduce by a standard continuation argument, since W(0) = and S 2 = o(e), that there exists an interval 
[0, jk], with To > fixed (independent of e, 5), on which 

W(t) <c 5 (e 2 + S 2 ) =o(e) 

which completes the proof of theorem [10] for the case of exact soliton initial data - part (i) of theorem [10] 



2.4.3 General initial data 

Part (ii) of theorem [TO] says that the behaviour described in part (i) also holds for nearby initial data: for a 
precise formulation it is necessary to consider the initial data for the gauge transform 
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Theorem 17 Under the same assumptions as theorem ] 1 (A let ^ be a solution to the Cauchy problem for \53\) 

hi 

, w, u, £) £ O sta b such that if we define x(x) = — (x — £) ■ a 5 (0, £), then 



with - ^ixt) € C(R;?i) and initial data ^(0) = (0(0), *0(O), A,,-(0), E., (0)) /lawng </ie following property. 



There exists A 



k 



-^•*(0)-^(0)-* sc , e (A) =o(e*). 

7~l 



(108) 



It follows that, if e is sufficiently small there exists Tq > 7 x(^ x ) anc ^ A(t) £ C 1 ([f J r3]j Ostab): a ^ as i- n 
theorem \l(A such that if^(t) is defined as in (|96p it satisfies all the conclusions of part (i) of theorem \l(A 



Proof It is only necessary to argue, as in the proof of lemma 1181 that under the stated conditions there 
exists A(0) e O s tab with |A(0) - A| = o(ei) such that W(0) = (0(O),<0(O), A^E^O)) = e" iexo • #(0) can 
be written as 

*(0) = (0sc ) e(A(O)) + 0(O) ! Vsc,e(A(O)) + ^(O),A i (O),E i (O)) , 

with 

(0(0), ^(0)) e N m 

where N\to) is the symplectic normal subspace, of codimension eight, defined in (|49[) . This is a simple 
consequence of the implicit function theorem, as is lemma 1181 There is only a slight modification required 



in that 0(0) 



0(0) depends on A(0), and so does ^(0), unlike the case considered in that lemma. 



However for small e this has no effect on the non-degeneracy condition required to apply the implicit function 
theorem. (Also the fact that xo grows linearly in x can easily be handled using the exponential decay in x 
of 4>sc,e,ipsc,e and their derivatives.) 

Now using | A(0) — A| = o(e^) we can deduce from (| 108[) that W(0) — o(e). Indeed for the electromagnetic 
components this is immediate since the gauge transformation leaves the electric field unchanged, and only 
shifts Aj by djXo, and this shift is put onto the background potential (and so does not contribute to W(0) 
since A is unchanged). The change of the electromagnetic components of the soliton induced by the change 
of A to A(0) are easily estimated in energy norm as 0(|A — A(0) | ) by lemmas 1551 and [Ml For the other 
components we just use phase invariance to estimate, e.g. 



0(0) - sc , e (A(O))|| 



L 2 = ||0(O) -e te *°0 SC) e(A(O))|U 2 
< ||0(O) -e ie *0 SC) e(A)|U 2 + ||e^ 
</so + 0(|A(0)-A|) = o(e*). 



0sc, e (A)-e ie *°0sc,e(A(O))|| i2 



From this point on, the argument can be completed as before: since (0(0), -0(0)) € -^a(o) is equivalent 
to the conditions (|99|) , theorems [14] and [16] can be applied to produce a local solution satisfying the growth 
estimate in i)2.4.2l 



3 Modulation theory 

In this section we state and prove some theorems which imply theorem 1141 which is needed in the proof 
of the main results (theorems [TO] and [TO)) . The proofs are a direct application of the developments in [21] , 
and so the presentation will be brief and reference made to [24] [16] for some of the calculations. The 
crucial point is that the conditions ([99]) are equivalent to a locally well-posed set of ordinary differential 
equations. Recall from ([47]) that, for e = 0, the soliton solutions are of the form (<fis,o, V'S,o)( x ; A) = 
e (/(»>(Z), (170; / W (Z) — 711 • Vz/o)(Z))J with A(i) an integral curve of the vector field Vq. Explicitly, the 
conditions ([M]) read 

(vM St0 (X)) L2 - (w,dT A <p s ,o(X)) l2 = (109) 

for A = -1,0, ...,6. 
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In the next two subsections we state two lemmas which prove that these constraints can be enforced 
thoroughout a time interval: 

• The first shows that by an appropriate choice of A(0), they can be assumed to hold in an open 
neighbourhood of the set of stable solitons in the phase space 7i. This shows that the class of initial 
data considered in part (ii) of theorem 1101 forms an open set containing the stable solitons. 

• The second shows that an appropriate choice of dtX implies that they are preserved for later times. 

3.1 Preparation of the initial data 

Lemma 18 Suppose that there exists X — (6,u>,u, £) € O s tab (so that (|39|) holds with u) = uj). Then, there 
exists e(A), n{X,e), such that, if \e\ < e(A) and 

= ||0(O) - <fec,e(A)|| H1 + ||V(0) - V^C,e(A)|| L 2 < K, (110) 

there exists A(0) € O s t a b depending differ entiably upon (0(0), "0(0)) such that 

(v(0),w(0)), determined by the first two equations of (|62|) at t = 0, satisfy the constraints (|109[) with 
A = A(0). Furthermore there exists C\ > such that 

|A(0) - A] + ||0(O) - 0sc, e (A(O))|| ffl + ||^(0) - ^?c, e (A(0))|| i2 < c{Hi. (Ill) 

Proof The condition in ([39]) allows this to be deduced from the implicit function theorem, see [24[ §2.3] or 
[IE] for details. □ 

3.2 Modulation equations and constraints 

Lemma 19 Let A(0) € O s t a b o,nd (v(0),w(0)) be as given in the conclusions of lemma [751 Let ^ be a 
solution to the Cauchy problem for \53\) on the time interval [0, T/ oc ] with regularity as in (|97[) . and such that 



sup 

[o,T ioc ; 



- * s ext (t) „ < N . (112) 



Fix a compact subset K of the stable parameter set O s t a b, which is the closure of an open neighbourhood 
of A(0). Then, there exists K2 > and T\ > such that, if \\{v(Q),w(Q))\\ H1( ^ L 2 < K2, there exists 
A(t)gC 1 ([0, T{\\ K) such that the constraints U09\) are satisfied for < t < T±, where v,w are as in ()62j) 
with ^ obtained from ^ via (I63p and (|96j) . The function t i— > A(t) is a solution of a system of ordinary 
differential equations (I116j) . 

Proof The proof of this is essentially the same as [241 § 2.5]. For clarity it is divided into three stages. 
3.2.1 Beginning of proof of lemma [191 

Equations (|66[) and ([57)) define a linear operator Ai\ in an obvious way: 

A4\(v, w) — (—dtv — ioj-fv + w, —dtw — iujjw — M\v) . (H3) 

and let A4* x be the formal L 2 (dxdi) adjoint of this operator. Then, by [2"4l, § 2.5], there exists an 8 x 8 matrix 
Dab such that 

Ml(-dZ^o, - D AB(-dT B ^s,o, dT B <l>s,o) + (IaJa)- (H4) 

B 

where the inhomogeneous terms V A are proportional to A — Vo(A): 

i A = i AB (x-Vo(x))B 
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with I AB smooth functions of x, which are exponentially decreasing as |x| — > oo; the precise formulae, 
which are unimportant here, can be found in 1241 §2.5]. A simple integration by parts then shows that the 
constraints in (|109p are satisfied on an interval containing the initial time, if they hold at that initial time 
and if the following is true 

(-dZ^s,o, h)v + {K A 4>Sfl,i2 +M) L 2 (115) 
- ihdx A ips,o, v) L 2 + (I A + ihdx A (t>Sfl, = 0, 

for all A = — 1, 0, . . . 6, and at each time in the interval. A calculation as in [24] . which is reviewed in the 
next stage of the proof in 8 )3.2.21 shows that these latter conditions are equivalent to the following system of 
differential equations 

(M(e)AB + 2ab(v, w, A)) (A - V (X)) B = F A (e, ¥ ext , A), (116) 

where M(e)AB is defined in (|117[) . 2ab is defined in (|118p , is given by (|122p and where the indices 
A,Be { — 1,0, 1, ...,6}, and we sum over the repeated index B. 



3.2.2 Explicit computation of the modulational equation (|116[) 

We write out explicitly the various terms in the conditions (|115l) . The first thing to note is that the overall 
expression is affinc in (A — Vq(A)) so we divide into the inertial terms, which are proportional to this quantity 
(and give rise to the left hand side of (|116p ). and the remaining force terms, which give rise to the right 
hand side of Q116P . The dominant contribution to the inertial terms arises from j[ , , while that to the force 
terms arises from j 1 ^ . 

To describe the inertial terms we need the following matrix, which, to highest order, describes the mass 
of the soliton: 

MUs(e) = (dxZipSfi, e^ B ^dx B 4>sc,e) l2 - (ft^s.o, e~ l ^dx B ^sc,e) ^ ■ (117) 
Then the dominant inertial term is 

(-^s,o,A)l» + (f£l<l>s,o,jZ)v> - M AB (e)(d t X - V (X)) B . 

Next, we have the following matrices, which may be thought of as corrections - owing to the presence of the 
perturbations v and w - to the "inertia" matrix above : 

Iab = (v, (J AB - idx B G c dx A i> s ,o)) - (w, (j% B + id XB Qcd\ A <f>s,a)) ^ ■ (118) 

We now present the abbreviations for the force terms appearing in the modulational equation. Firstly, 
we have what is effectively the Lorentz force term. 

F L A = (fr A il>s,o,ff)v> - (9xI<fe,o,iV 

'dx^sfiiiea^f^^^ ( 119 ) 

'd\ A <t>sja, iea ' x (ij(u - ea w , e ) - u • V) f u>e + eR/ u>e 
We also have a force + F^ due to the nonlinear interactions, where 



= -{d XA cbs.oM) L2 , (120) 
K = (&s,o,Ji+^( 7 « w , e + a^ + i ) W ) i2 (121) 
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We abbreviate the total force as follows: 

P^=P^+P1 + P^. (122) 
Bound for the inertia matrix. It follows from the definition of Iab that 

\1ab\ =0 (!!(«,«;, A, E)||«) . (123) 

Bounds for the forces. Firstly, the main force term can be bounded as 

Fi = 0(e), (124) 

because of (|187p . (|188p and (|154[) . For some values of A there are better bounds: 

K = 0(e 3 ). (125) 

Referring to (I119f> , and using lemmas IA.1.31 and [3JJ we deduce that 

Fq = 0((?)-(if u ,eRf u ) L * 

+ (detps,o,iea S ,x f UJ ) L 2 - {dg(j)s,o,iea^ x {iuj^f u - u • V/ W )) L 2. 

By the reality of / w and the Coulomb condition the last three terms vanish, proving the bound (|125[) . Also, 
for A = 3 + j we have an improvement: 

Ft +1 =0(e 2 + e5). (126) 
To establish this, we first argue as above that 

Ff +J - 0(e 3 )-(cWs,o,eRX) i2 

+ (d u jips,o,iea ' x f LJ ) L 2 - {d u , <f> s ,o, iea ' x (iuj-f f u - u- V/ W )) L 2. 

Now referring to the formulae in lA.1.4l we see that d u j(j>s,o = even + iodd, while d u jips.o = odd + ieven where 
even (resp. odd) means a real valued function which is even (resp. odd ) as a function of Z. The bound 
asserted then follows by inspection and use of lemma [37] 

Next, (JSS} implies, by ([1501) ), (BSD > and by the Holder and Sobolev inequalities, that 

|F2| = O (e 2 \\(v lW , A,E)||„) +0 (\\(v,w, A,B)\\h + \\(v,w, A,E)||^) . 

Finally 

|FP | = O (e|| («, w, A, E)||„ + e||(u, to, A, E)||^ + e 2 ||(i;, to, A, E)||^) . (127) 

This is obtained directly from the formula above by means of the Sobolev and Holder inequalities and using 
the bounds in qA.2.1l and 3A.2.21 

3.2.3 Completion of proof of lemma [T9l 

The matrix M(e),4B is invertible for small e on account of the stability condition ([55| and lemma ([35]). Also 
the matrix Iab is small when (v, w) is small, so that in this case the system of evolution equations (|116[) can 
be manipulated - as in the proof of theorem 2.6 in [53] - to form a system of equations of the form 

X = V (X)+V i (e,¥ ext ,if,X). 

This is almost a locally well-posed system of ordinary differential equations - there is a slight modification 
of the standard proof from [24] required: VP is known to exist already, but (v,w), determined as in the 
statement, depend on X(t) through the gauge transformation (|63[) . which is nonlocal in the £ component of 
A, and so V\ is similarly nonlocal. To allow for this it is necessary to augment A by the nonlocal quantity 
appearing in ([6"3"]) . which is in fact x(t, £). Call A = (A, x(i, £(*)))) then there is a locally well-posed system 
of ordinary differential equations of the form A = W(A), e, ^ s ext , allowing the proof of lemma IT51 to be 
completed in the same way in [24]. □ 
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3.3 A bound for A 

Lemma 20 In the situation of the previous lemma, 

|A- F (A)| = O (e+ \\(v,w, A,E)|& + e\\(v,w, A,E)||«) 
in the limit of e going to zero. 

Proof The function A(t) is obtained as a solution of the modulation equations (|116p . Referring to the 
bounds for the inertial matrix and forces in H3.2.21 it is immediate that for e, |j (v, w, A, E)||^ sufficiently 
small the bound claimed holds. □ 



4 The Lorentz force law: proof of theorem [12 

The starting point is (|116|) . Define 

M AB (0) = (d^ips,o, d\ B 4>Sfl) L2 - (dT A ^s,o, d^^s.o) , (128) 

and observe that by lemmas [35] and [34] M.ab(s) — Mab(0) = 0(e 2 ). Using this, and referring to the 
decomposition of Fa in equation (I122p . and the associated bounds following it, we infer that 

(m(O)ab + 0(e 2 + W^(\-V a ) B = F L A + 0(eW? +W), (129) 

where is as in (|119p . Since the right hand side is known, up to the stated error term, it is now just 
a matter of calculation to obtain explicit forms for the left hand side of these equations, and thence to 
deduce theorem Q21 The calculation is done in [2H §A.7], using a set of functions defined in ^ A. 1.41 which 
are convenient linear combinations of the dx A ((ps,o,ips,o)- We now record the conclusions. 
Using (|102p . the A = component of (|129|) reads: 

WfJlip = F^ + 0(e 2 ) + 0(W), 

with a formula for Fq given in {TTal) which indicates that F^ = 0(e 3 ) (see §3.2.2jl . and all together: 

dUu\\fu\\ 2 L 2)u = 0(e 2 ) + 0(W). (130) 

Similarly, the bound (|126[) for Fg +J implies the following equation for the centre of the soliton: 

£ = u + 0{e 2 ) +0(W) + 0(eS). (131) 
Next, using (fT3Q)) and (fTU^j) . the A = i e {1, 2, 3} component of (fT2"9")) reads 

di 

again with Ff given in (| 1 19[) as: 

Pf = (^s,o,iea 5 ' x f t 

- (^d^4> S:0 ,iea S ' x - ecw) - u ■ V) f^ e + eR/ w , e ^ , 

where the operator R is defined in (fHTj) . Here, on the left hand side, ||/ w ||^2 = / /w(Z) 2 e? 3 Z and by the 
Lorentz transformation (|4"2"]) d 3 Z — 7d 3 x. The inner products on the right hand side are in L 2 (d 3 x). It 
remains to simplify this expression for Ff : firstly, 

(%V'S,o,«eao X /",e) i2 ~ (d^4>sfi,iea^ x {i-f(uj - ea^^e) -u-V)/ w , e ^ 

= (%V's,o^eao X / w ^ 2 - ^%0s,o,ieao' x (i7w-u- ^)fJ) 2 + 0{e 3 ) 



[ '\ l|V/.|| 2 L2 + lo 2 \\fJ 2 L ,h^\ = Ff + 0{W) + 0(e 2 ), (132) 



L ' 2 (133) 
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by lemma 1551 



= ((i 7 w - u • V)/ w , ieVaJ*/,,,) + 0(e 3 ) 

by integration by parts, 



- eu\\f u \\l a [V ia d (t, £) - k d (t, |) - u • Va d (i, £)J + 0(e<5 + e 3 ), 
by JBSJ and lemma[351 (Again, ||/<j£ a = / / w (Z) 2 d 3 Z.) But also, referring to ([51]), 

(-%^,o,eR/ w , e ) i2 = 7 we J f 2 (Z)Vu ■ a 6 (t,x)dx, 

= we|l/«lli»UiV r af(t,€) + 0(e*), 
again using lemma l38l Adding together these contributions, we end up with 

F L = ew||/ w ||£ 2 (Vag - (9 t a 5 ) + u x (V x a s ))(t,£) + 0(e 3 + eS), 

which is the required form of the Lorentz force law, as given in theorem [TSl once we note that J ojf 2 = 
j{u-e a )fl e + 0{e?). □ 

5 Proof of the main growth estimate 

In this section we are concerned with the proof of theorem 1 161 In order to control W it is helpful to introduce 
a quantity W which allows us to take advantage of certain cancellations occuring in the energy identity to 
handle some of the nonlinear interaction terms which would otherwise be difficult to estimate directly. The 
direct nonlinear interactions between v and A arise from terms in the Hamiltonian obtained by expanding the 
expression | J | (V — ieA)<p\ 2 in terms of v, A by means of (There are also indirect interactions mediated 
by Ao via the Gauss law, but these are easier to estimate.) In H2.1.3I this expansion of | J |(V — ieh)<f)\ 2 
is carried out explicitly, and, including also the quadratic part of the Taylor expansion of the potential V, 
leads to the introduction of the quantity: 



H(v,A) = Y,H {n 

n=2 

= -Je(( w ' A )' b( " _1) ( v ' A 



n=2 



where the superscript n (resp. n — 1) indicates the homogeneity in v, A of the term in the expanded 
Hamiltonian (resp. of the term B^™ -1 ) in the expanded evolution equations (|67p . (p9")0 ; see ^2. 1.31 for explicit 
expressions and explanations. Using these definitions we have an alternative form for the expanded evolution: 
equations (|66|) . (|68p can be written in the form 



d t (v,A) = (w,B) - (i(juj + h)v,0) - (d t X — Vo(A)) • {d x <t>sc,z,dxA sc ,e) + + (*n,0), (134) 

with $11 — j^ 1 . The remaining two equations (167p . (|69p can be written: 

d t (w,E) = (-D V H,-D A H) - (i(jLj + h)w,0) 

- (d t \-Vo(X)) ■ (lhi>sc,e,dxEsc,e) + (j§,Q) + ($21, $22), (135) 
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where h is defined in (j72"|) , and $21 = Jjj + A/", and $22 = j| 7 are given in terms of the inhomogeneous terms 
defined in ^2 . 1 .31 notice that the inhomogeneous terms 32 11 1 J 1 ^ > 3& 11 > J 1 ^ are included in the first term on 
the right hand side of (|135|) . 

To study these equations it will turn out that the following quantity is useful: 

W= ^\\w- itHIl* " ^^IMlL + + ((«>,% u- V(u,A)^ +H{v,A). 

We can think of W as follows: it is formed by adding to the Hessian of the augmented Hamiltonian W 
those terms arising in the expanded Hamiltonian (when we input the perturbed solution ansatz (|62| ) which 
describe the interactions of the fields (v, A) with themselves and with the external electromagnetic field. An 
important reason for introducing W is that the following two lemmas imply a long time bound for W, and 
hence a stability estimate in energy norm. 

Lemma 21 In the situation of theorem \15[ so that 

• A lies in a compact subset, K C O s t a b, 

• (v,w) satisfy the constraints (|99p . and 

• W is equivalent (uniformly onH) to \\(v, W, A, E)||L, 

assume that W < 1, and that e = o(l) and e — 0(6). Then, there exists a constant c(K) > such that, for 
all A e K, 

cW < W < -W. (136) 
c 

Proof Referring to the formulae in pXBI for the which occur in the definition of H, it is a straghtfor- 
ward consequence of the Holder inequality that 

2 2 

W = W + 0{- c W) + 0(%W) + 0(^-Wh + 0(e 2 W) + OieWh + 0(e 2 W 2 ), 
o z 

lemma [57] and the assumptions on the external field in HI. 21 The lemma follows immediately. □ 

Notation 22 In the following we write, f = 4t (O(A) + o(B)) if there exist C 1 functions g, h such that 
f = d_( g + h) and g = 0{A) and h = o(B). 

Lemma 23 Assume the hypotheses of theorem \lb\ It follows that, 



dt 



d ■ ■ • ' • -' 



= Jt 1° [ eWi ) + "WJ + ( e4 + ( e + j) W + ( e2 + eS ) Wi ) . 1 1 : > 7 ' 



in the limit of e and W going to zero. 

Proof See SJEU □ 
5.1 Proof of theorem 1161 assuming lemma [53] 

Proof Integrating up equation (|137[) . and using the Cauchy-Schwarz inequality, lebW 1 ! 2 < +eS 2 + eW, we 
infer the existence of a constant c > such that, for t £ [0, T2/e], 

,2 , c2 



W(t) - W(0) < c I e 2 +S Z + |e| J W(s)ds \ , (138) 

as long as e = O(S). By Gronwall's inequality and lemma [2"31 for \e\ sufficiently small there exists a constant 
c > such that, on [0,T 2 /e], 

W(t) < c (W(0) + e 2 + 5 2 ^ exp[c|e|t]. (139) 
By lemma [2T1 the result is proved. □ 
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5.2 Proof of lemma 

5.2.1 Beginning of proof of lemma 

By the assumptions of theorem [TCI we have a solution of equations (|134[) . (|135p satisfying the conclusions of 
theorems IT4l and [TCI so that the constraints (|109p hold and W = 0(e). Then, by lemma I2T1 and theorem [TCI 
there exists c > such that 

-W < || (v, iu, A,E)|& < cW. 

Also since W = 0(e) the bound (|102jl holds, and will be used in the course of the proof. The estimate for 
W will be obtained as a consequence of the energy identity for (|134|) . (|135|) . so the next stage is to write that 
identity down and separate the terms out in a way that allows them to be usefully estimated. 

5.2.2 The energy identity for f[66 jl - f|69 ]) 

^-W = (d t (w,E),(w,B)+u-V(v,A)-i'yw(v,0)) (140) 



d t (v,A),(-D v H,-D A H) -u- V(ui,E) +iju(w,0) 



L 2 



+ / d t h(v, A)dx - dt(ju) (iv,w) L 2 + (d t u- V(u, A), (io,E)^ 



L 2 

Here we have introduced a notation h(v, A) for the integrand defining H , i.e. 

H(v,A) = [ h(v,A)dx (141) 



~ [^((vM^H^A.)) dx. 

J o 



Explicit expressions for the nonlinear operators !>(™ -1 )(i>, A) show that they depend on t, x, and the dji in 
the final line of (|140[) refers to differentiation with (v,A) held fixed; similar conventions will be understood 
below. 

Substituting for the time derivatives from (| 134|) and (| 1 35[) , and noting the usual cancellations which occur 
in the derivation of the energy identity, we obtain the following expression: 

—W = Q1 + Q2 + Q3- (D v H,ihv) L 2 + J (d t + u • V)h(v, A) dx (142) 

- (iv, (u • Vh)w) L 2 ~ dt^uj) (iv,w) L 2 + (d t u ■ V(v, A), (w,E)^ 



L- 



where 



Qi = (($21, $22), (w,E) + u-V(u, A) -i7w(«,0) 



L 2 



^u.O),^,-^) + u. V(w,E) - i7w( M) 0) 



L 2 



\ / L 2 

- ( (j?, Js) . (~ A,#, ~ d aH) ~ u • V (w, E) + (ib, 0) 

\ / L 2 

and Q 3 = -(<9 t A - Vb(A)) • Q3, where 

Q 3 = ({K^Psce, dxE sc ,e) , (w, E) + u • V(u, A) - i 7 w(t; 1 0) ) 

' (Kfcce, dxAsce), {-D V H 7 -D A H) - u • V(u>, E) + z 7 ^(w, 0) 



L 2 
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We control Qi, Q2, Q3 in the next three subsections before completing the proof of lemma |2"B1 In the course 
of estimating the various terms we will use bounds for TV, h and the 4>'s (which may be read off from those 
in 92.1.40 . and the bounds for A Q in 9A.2.21 



5.2.3 Estimation of Qi 

The following proposition is the main result about Q\ needed for the basic growth estimate: 
Proposition 24 In the situation of lemma[£3\ 

Qi = d t (o{W) + 0(eW^fj + o(e 4 + eW + e 2 W? + eSW^ . 
Proof Substituting from (] 1 34[) and (| 135[) we obtain: 

Q 1 = {d t \-Vo)-(d x Asc,e$22) L 2 (143) 
+{d t X - Vq) ■ {(dx<i>sc,e, $ 2i) L2 - (dxi>sc,e, $11, ) 2 
+ (ieA (ij(uj - ea^.e) - u • V)/ W , e , $11) 

-(Mo/c, e ,Vi ), ($21^22))^ 

+ ((d t + u • V)v, $2i) L 2 + (ihv, $2i) L 2 - (ihw, $n) L 2 
+((d t + u • V)A, $ 22 ) L2 - ((d t + u • V)w, $n) i2 

since $12 = 0. 

Estimation of the first line in Q\ The first line of Q\ is easily seen to be small, since $22 = —e 2 ^ S ' x fui.e is 
0(e 2 ) in every L p by the bounds in 92.1.41 Together with the fact that, IISaAscbIIlp = 0(e) for p > 3, by 
([STj) and the results of appendix IA. 1.21 this implies that (d\Asc.e, $22) l 2 = 0(e 3 ), and so by (|102[) the first 
line is 0(e 4 ). 

Estimation of the second line in Q\. The second line is smaller than appears due to a cancellation which is 
a consequence of the modulation equations, pi5[) or pi6[) . To see this, we refer to the decomposition of the 
force on the right hand side of (|1 given in 93.2.21 and using the definitions of the $7 j in (|134[) , (|135[) , we 
see that 

'^so, $21} l2 - (^I^sc, e ,*ii) i . a = -F5-FS + £rr A 

= - (M(e) AB + Iab) (A - M)) B + + Err A 



where 



Err A = (d\ A (j)sc,e - d\ A ^ S fi, ^21 j - (d\ A ip S c,e - d\ A ip s ,o, $11 f ., 



Using lemma EH the bound ((95]) for Af, and the fact that from 92.1.41 <5>n = j{* and $21 - AT = $ arc 
0(e), we deduce that |£Tr\4| < ce 2 (e + IU + 1-U 5 / 2 + e 2 W x l 2 ). Next notice that lemma 1351 implies that 
M(e)AB -M(0)ab = 0(e 2 ). Therefore since M(0)ab = -M(0) B A the largest term drops out and the second 
line of Qi can be rewritten as 

(M(e) AB - M(0) AB + Iab) (A - U ) A (A - Vb) B - (F* + Err A ) (A - U ) A 
which, by the above and (fll^ . (fmf is 0(e 4 + e 2 ^ 1 / 2 ), for small e and VU. 

Estimation of the third and fourth lines in Qi. Using lemmal39" l (|95 j) . the bounds in 92. 1.41 and the properties of 
in appendix lA.1.11 the third and fourth lines can be estimated immediately to be O (^W 1 ^ 2 + e 2 W 3 ^ 2 ^j . 
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Estimation of the fifth and sixth line in Q\. This requires care because h is unbounded as a function of x. 
This makes it essential to separate the nonlinear term Af in $21 from the other terms (which are exponentially 
decreasing in x and can thus absorb the unboundedness of h) . Therefore we estimate first of all the quantity 

(ihv, $21 - M(f u , ei f u , v)) L 2 - (ihw, $n) L2 - O (e 2 W^ , (144) 

by (|102p and the bounds for h recorded in i)2.1.41 Next, write the first term on line five, together with the 
missing piece (ihVjAf)^ from the previous estimation, as the sum of two quantities: 

{(dt + ih + u- V)v,N) L 2 + Rem, 

where Rem — ((d t + u • V) v, $21 — N)l 2 - ^ ' s sri0wn m lemma |4"D1 that the first of these quantities is 
dt(o(W)) + 0(eW + e 3 W^). To complete the proof of proposition [24] we need to estimate the sixth line 
and the quantity Rem defined above. This is done by means of the integration by parts identity (1195|) . and 
taking advantage of the fact that 

(d t + u- V)/ W>e = (A — Vb(A)) ■ d x f UlB , (145) 

is 0(e) by flTU2J). Together with fTM]l . this implies that 

\\(d t + u- V)<t> u \\ LP = 0(e(e + 5)) (146) 

for all p and all I J except for I J = 21; but in that case (|146[) holds instead for $21 — M = j| i \ (which is 
what is actually needed to estimate Rem). Putting this information into (|195|) . we infer that the sixth line 
and Rem are dt(0(eW^)) + 0(e(e + 8)Wi), which is sufficient to complete the proof of the proposition. □ 

5.2.4 Estimation of Q 2 

The terms in Q2 arising from j®, can be estimated in a straightforward way by the Holder and Sobolev 
inequalities, because of the exponential decay of /^. e , and using lemma [3"51 to bound Aq. For example, 

(^w - ijujv + u • Vv, ieA (17 (w - e« u , e ) - u.V) fu,ej 2 = O (e 2 W) (147) 

by Holder's inequality, since f u , e and V ' f u , e ar< 3 bounded in every LP norm and ||j4o||lp = O (eW^j for 
3 < p < 00. For the terms involving j° = 'S/Aq we can estimate, 

(u ■ VE, Vio) = (divE, u • Vi ) i2 = O (e 2 W^ , (148) 

since ||divE|| i3 /a = O (eW^j and ||VA ||l 3 = O (eWi^j . Consider next the terms ? . Refer- 

ring to the explicit expressions for D^H given in ^2.1. 3[ starting with ([88]) . we see that the resulting terms 
can all be estimated in a straightforward way (using the bounds for S7A in appendix lA.2.2[) to be 0(e 2 W), 
except for one, namely: 

<AA,Vi ) L 2, 

but this vanishes by the Coulomb condition, and so Q2 — 0(e 2 W). 

5.2.5 Estimation of Q 3 

The quantity Q3 is smaller than it appears due to the constraints. To see this first recall that, as used 
above already, ||5AAsc,e||iP = 0(e) for p > 3, and ||9AEsc,e||i> = 0(e) for p > 3/2, by (|5Tj) and the 
results of appendices IA.1.21 and I A. 1.31 Referring to the expressions for D^H in ^ 32 .1.31 this means that the 
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electromagnetic contributions to Q3 can be bounded as 0(eW? ). But also, the expressions for D V H in i j2.1.3l 
imply that 



dxcbsce, —D V H + M x v) l2 = O [eW* 

Therefore, up to O(eWi), we deduce that Q3 is equal to 

(u ■ Vw - iujjw - M\v, d\(f>sc,e, ^ - ■ Vw - iuJjv + to, d\ipsc,e} ■ 

Now the identities in appendix IA.1.4I and the constraints (|109[) imply that this expression vanishes if 
( l > sc,e,'4>sc,e are replaced by 0s,Oi' ! /'S,o- But by lemma l3"5l this can be done at the expense of an O fe 2 W^ 

error. Therefore, since (A — V ) — 0(e) by (|102p , we deduce that Q 3 — O (^e 2 W^ . 
5.2.6 Completion of proof of lemma l23l 

The previous subsections have provided the requisite information on the Q's, and so it now suffices to 
control the remaining quantities in (|142p appearing after the Q's. The following two propositions treat the 
two quantities on the first line of (|142|) . 

Proposition 25 Assume the hypotheses of lemma \23\ . It follows that, 

J{d t + u-V)h(v,A)dx = -iy^(>,,A),(d t + u- W)M^(v,A)) dx 



n=2 
() [r]V+ 



Proof Observe 



the fact that a S ' x is pointwise O(y), but its derivatives are O(l), in particular ||(9 t + u • V)a d,x ||ioo < 
2(||a|| L - + ||Va|| £ -). 



• the identity 

(Ot + U- V)f u , e = (A — V ) ■ d X f u , e , 

which shows that the left hand side is 0(e) in every L p , by (|102p and the exponential decay properties 
in appendix IA.1I Similarly, \\(dt + u • S7)a UJte \\w 1 ^ is 0(e 2 ) by (|102[) and the bounds for a Wje in 
appendix IA.1.21 

To prove the proposition now, just use these observations to estimate with Holder's inequality each of the 
terms arising from differentiation of the expressions for in i)2.1.3l □ 

Proposition 26 Assume the hypotheses of lemma \2S\ . It follows that 



(D V H, ihv) L 2=0 l^eW + —W 

Proof Using the notation in (|88[) for the Frechet derivative D V H, we have 

\(D v H,ihv) L *\ = \(B(v, A), (ihv, 0)) L 2\ (149) 

and we can estimate term by term, but some care is needed since h is unbounded as a function of x, see 
(|72[) . In addition to the first point in the proof of the previous proposition, we use the bounds for h recorded 
in i)2.1.41 Those terms in (|149[) arising from lK 3 ) vanish identically, while of those arising from B^ 2 ^ the 
only non-zero ones are proportional to e(hv, VvA)^ . By the Coulomb condition and the bound for V/i 
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from EH this term is 0(e 2 W 3 / 2 ). It remains to bound those terms arising from '. Of these, it is 
straightforward to bound those arising from B 12 as 0(eW) by the second fact just mentioned, and the same 
goes for those arising from M\ in In = — M\ + eR + S. However, there is a single non-zero term arising 
from eRu which is proportional to 

(hv, a s ' x ■ Vv) 

2 

which, with an integration by parts, can be bounded as O(^-W), but, again, only after taking into account 
the Coulomb condition V ■ a S ' x — 0. Finally for the terms arising from S we see from (|5Tj) that 

(ihv, Sv) L 2 = ej J [2a (i)ie (wu • Vu) + u • Va^^lwl 2 ] dx = 0, 

so that {ihv, Sv)l2 = 0, and the proof of the proposition is completed. □ 
The remaining terms on the second line of formula (|142p are easily estimated as 0(eW) by (|102|) . and 
the proof of lemma [23] is completed. 

A Appendices 

A.l Further properties of the solitons 

A. 1.1 Exponential decay properties of the solitons 

The e — solitons in the nonlinear Klein-Gordon equation (|23p are exponentially localized: to be precise we 
have the following estimates for the profiles functions fu,9uj'- 

lim sup V V a f u Exp[\x\ (V«i 2 -oJ 2 - e)} < oo Ve e (0, Vm 2 -Lu 2 ), (150) 



\x\ — »oo 

l«|<3 



together with 



and 



lim sup ^/ a g u Exp[\x\ ( V m 2 - lo 2 - e)] < oo Ve € (0, \/ m 2 - w 2 ), (151) 



\x\ — >oo 

|q|<3 



f 



lim if- = -^m 2 -UJ 2 , (152) 

while Ve > 0, there exists c(e) > such that 

f u (\x\) > c(e)Exp[-\x\ [^/m 2 - lo 2 + e)]. (153) 

(See theorem 1.4 in [24]). Exponential decay also holds for the solitons coupled to electromagnetism for 
small e: 

Lemma 27 Suppose that |e| < e±, for some e\ > 0. Under conditions &29W3ty) on U , 

\D a f u , e (x)\ < CEx P [-k\x\] (154) 

for positive constants C and K, and where a is any multi-index with \a\ < 2. Furthermore, the constants C 
and k are independent of the coupling constant e. 

Proof See [17]. □ 
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A. 1.2 Some estimates of the soliton electromagnetic potential a 
Lemma 28 For each f G H 2 (R 3 ), there exists a unique a € iJ r (R 3 ) such that 

- Aa + e 2 fa = ujef. (155) 

Furthermore, the map A : iJ 2 (R 3 ) — ► H (R 3 ) defined by A(f) — a is continuously Frechet- differ entiable. 
Proof This follows from standard arguments. □ 
Lemma 29 Suppose that f G H 1 (M 3 ). Suppose further that a solves 

- Aa + e 2 f 2 a = eujf 2 . (156) 

It follows that Va^Wa G L 2 (R 3 ) for any i,j G (1,2,3). Furthermore, ||V i V- 7 'a|| i3 , ||Va|| L2 , \\a\\ Lco = 
0(e) 

Proof 

r \Va\ 2 +e 2 f 2 a 2 = ecu I f 2 a (157) 



from which it easily follows via Sobolev's inequality that 

||Va|| L2 <ce ll/H^II/1^3. (158) 
Next, since -Aa = e (u) — ea) f 2 , we have 

|| Aa|| L2 < e (w \\f\\l< + e \\a„, e \\ Le ||/||* B ) . (159) 

By the Calderon-Zygmund inequality, we have that for any i,j G (1,2,3), 

||V l V J a|| L2 = 0(e). (160) 

By Sobolev's inequality, we have thus shown that a G W 1,6 and hence by Morrey's inequality, ||a|| L oo = 0(e). 
□ 



Corollary 30 Suppose that / W)B G iJ 2 (R 3 ) solves 

- A/^e + m 2 fiu,e ~ (u>- ea Ui e) 2 /w,e = 0(fu,e)fw,e, (161) 

where a Ul e G i? r (R 3 ) is a non-local function of f u , e uniquely determined by 

- Att^e + e 2 /^ e a: W!e = we/ 2 e . (162) 

TTien, /„. G tf 4 ( 



Proof Differentiate the equation for f ue and apply the Calderon-Zygmund inequality. □ 
This leads naturally to the following lemma. 

Lemma 31 Suppose that f G if 4 (R 3 ) and that a solves 

- Aa + e 2 f 2 a = eujf 2 . (163) 

It follows that Va G W^ P (R 3 ) for any p G (§,oo). 

Proof Differentiate (|163[) . and apply the Calderon-Zygmund inequality (using the Holder and Sobolev 
inequalities if necessary) to get the result. □ 
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Lemma 32 Suppose that f £ H (K J ) and that a solves 



Aa + e 2 f 2 a = euf 2 



(164) 



It follows that 



where sgn(x) — x/\x\ for x ^ and sgn(0) = 0. 



(7) 





UJ 


) a < 






e 



Proof Assume that / in C^°(R 3 ). Define a + = max(a,0) and a~ = max(— a, 0). Suppose cue > 0, then 
by a weak maximum principle (theorem 8.1 in [8]), a > 0. Now, .Ag = a— — solves — AA + e 2 |/| 2 A = 0, 
therefore yl < by the same weak maximum principle. Hence, < a < — . Similarly, if —we > 0, then 
> a > — — so that ||a|| ioo < |^|. The lemma follows by approximation. □ 

Lemma 33 Suppose that f^^ e and a w , e are as given in theorem^ Then, 

; da u 



V 4 V J - 



0(e) 



LP 



(165) 



for p G (1, 00), and i, j = 1,2,3. In addition, 



dX 



0(e) for any p £ (|, 00) 



Proof From lemma l28l and theorem [6l is a well-defined object. We note that 



A 



dX A + ' J - f ' ' ^7a7 = e/ <^-^ + 2e/ "' e(w " ea - e) "^ 



(166) 



from which 



A 



d\ A 



L*' 



0(e) for p £ (l,oo) follows immediately. The lemma follows trivially from 
repeated differentiation, the Calderon-Zygmund inequality and the Holder and Sobolev inequalities. □ 
Let ((x;\) be the unique solution in H 1 of (fS"2l . —AC = — 711 ■ \/a u , e (Z), which takes the Lorentz 
transformed solitons into Coulomb gauge. Then 

Lemma 34 



IIW'cl 
|v*v J *s A c| 



Lp 



LP 



0(e), 
0(e), 



(167) 



for p £ (|, 00) and i, j — 1, 2, 3. 

Proof By (|52p . and its derivative: and 



a d > „ d 

-A— — ( = -7U • V-T— a w ,e 



the result follows by means of lemmas [35] an d [ 



□ 



A. 1.3 Differentiability 

Lemma 35 Let f u>e £ H 2 be given by theorem^ Then it is a differentiate function of uj and satisfies, for 
small e: 

||/ w> e - U\\ H 2 + \\8 u f u>e - d^fJ H2 = 0(e 2 ). (168) 
Proof See El . □ 
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Lemma 36 Let h u = h w — uq Ul where h w = ff($g, e (0,w,0,0)) while q^ = Q($s, e (0, w, 0, 0)). Then 

-j-h u = —q u . (169) 
aid 

Proof Following the argument given in j!0| , we note that 

■jJL = -fe + (h'($sA*o)) ~ ^Q'(%,e(A )), ^ $ s,e(A )^ ^ , (170) 
where A = (co, 0, 0, 0). The result follows from the fact that H'($> s>t ,(\ )) - wQ'($s :e (A )) = 0. □ 



A. 1.4 Some identities involving (^dx4>s,o, d\ipS,o 



The explicit calculation of the modulation equations can be carried out by making use of the following 
functions (a v i,(Z(x ) A); A), 6^(Z(x, A); A) from [24] : 

6_i(Z;A) - .g L ,-^u.Z^, (171) 

6 (Z;A) = */„,, (172) 

h(Z;X) = V l z / w (Z), (173) 

6 3+4 (Z;A) = i V J z ^(Z)-z^7((7P u + Q u )Z)J w (Z), (174) 

while 

a_i(Z;A) = -7" 1 6 + (7U.Vz-i7^)^-i, (175) 

a (Z;A) = ( 7 u.V z - tjuj) b a (176) 

a;(Z;A) = ( 7 u.V z -«7w)6 l , (177) 

a 3+l (Z;A) = (7P U + Qu) Z)«6j + (tu-Vz - *7") &s+i, ( 178 ) 
where i, j = 1, 2, 3, 5^ = and 

Cm - 7 2 (u • z)(p u ) ri + Ij=jl( u • z)(Q u ) Jl + 2_i(Q uZ ) 4U3 , 

These are convenient for computation of the modulation equations because the linear span of the d\ A (<f>s,a, fps,o) 

is the same as the linear span of the (bA, — oa)-(To be precise: except for A = j e {1,2,3}, we have 

dx^sfi^Sfl) = (b A ,-a A ), and for A = j we have {4>s,o, ips,o) = ~(lPu+Qu)jk {h, -a k )+^u j (b , ad).) 

The following identities are equivalent to lemma 2.2 in [24], and can be obtained by differentiating the 
Euler-Lagrange equation Fq — 0, where Fo is the augmented Hamiltonian ([4"5|) : 

(iju- u- V)9 Ao 0s,o - d Xo tps,o = 0, (179) 

(i7W-u-V)&s,o-M A ^s,o = 0, (180) 

(i-fuj - u ■ V) dxj^sfi ~ dx 3 ips,o = 0, (181) 

(ijcj-u- V)^s,o- Mx^s,o = 0, (182) 

(«7W - u • V) 9a_! (t>s,o-dx_i ips,o = — d Xo <ps,o, (183) 

7 

(«7^ - u • V) dx^ipsfi - M x dx_ 1 4's,o = — d Xo ips,o, (184) 

7 

(i7w - u • V) dx 3+j <ps,o - dx 3+] ips,o = -9x^3,0 - 7uUjdx <ps,o, (185) 

(17^ - u • V) dx 3+J ips,o - M x dx 3+j 4>sfi = -dxjipsfl - iuujdx ips,o (186) 
where the index j runs from 1 to 3. 
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A. 2 Some estimates 

A. 2.1 Estimates related to the external field 

Lemma 37 Let f be a measurable function with (1 + \x\)f S L 1 . Then if a S ' x is as in (|64p 

e4 x fj Lf <ceL 1 \\(l + \x-^\)f\\ LP , 
and 

||ea^/|| iP < ceii ||(1 + \x-£\)f\\ LP 
for p S [1, oo] . If in addition f even is an even function of (x — £) and (1 + \x\) 2 f even £ L 1 then 



with L\,L 2 as in |7[). 

Proof Recall HMD and Writing 



< cL 2 S II (l + \x-£\) 2 f, 



evert 1 1 j^i 



(187) 
(188) 

(189) 



og(t, x) - oS(i, e) = (x - • / Vag(*, £ + fl (a; - £))ds 



(190) 



etc, by the fundamental theorem of calculus, the result then follows, using the fact that the gradients of 
a^a" 5 are bounded independent of 8 by assumption (see £)1.2[) . For the proof of (|189[) . it suffices to use 
the identity for Va^ corresponding to (|190[) . and then substitute this back into (| 1 90[) and use the fact that 

/(X - £)/««, = 0. □ 

Similarly, we have the following bounds: 



Lemma 38 



Ki + is-eirHft+u-vja''* 



< d(l*l + H) 



(191) 
(192) 



w/iere we use pM]l . ©, and Ci = Ci(ii,i 2 ) and C 2 = C 2 {L 2 , ||(1 + |x|)/|| L i). 

A. 2. 2 Estimates for the time component of the electromagnetic potential 

Lemma 39 Given (v,w) £ H 1 x L 2 and A G O there exists a unique A$ £ H 1 solving (|7ip such that 



for p £ (1,3]. Consequently 



VA 
A 



L" 



O ( e|| («, to, A, E) ||« + e|| («, A, E)||^ 



(193) 



L'J 



satisfies the same bound for 3 < q < oo 6w Sobolev's inequality. 



Proof From Gauss's law (|7"Tj) . we have explicitly 

- AA = e (i/ w , e , w) + e (iv, (47(0; - ea w , e ) - u.V) + w) 
By Sobolev's and Holder's respective inequalities, 



(194) 



AA 



L" 



0(e\\(v,w,A,E)\\ H + e\\(v,w,A,-E)\\^ 



for q G [1, |]. The lemma follows from the Sobolev inequality and from the Calderon-Zygmund inequality, 
HO section 9.4]. □ 
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A. 2. 3 Integration by parts and simple averaging 



First we recall the phenomenon of averaging in the context of ordinary differential equations, in the simplest 
possible case of the perturbed harmonic oscillator. Let g be a C 1 function of t 6 K, with \g\ < M and 
\g\ < N. For < e << 1 let y e be the solution oiy + y = eg(et) with initial data y e (0) = yo, y e (0) = y\ (fixed 
independent of e). Then y e — y° is O(e) in C 1 ([— T, T]) norm for times of T = 0(±). One way to prove this 
is to define / = e^ 1 {y e — y°), which solves / + / = g(et) with zero initial data. Let E(t) = (f 2 + / 2 )/2 be 
the energy; it satisfies E(0) = and E(t) = f(t)g(et). Now an integration by parts gives 



which, by Gronwall's inequality, implies E(t) — 0(1) for t = 0(-) as claimed. To conclude, this simple fact 
- that a small slowly varying inhomogeneous eg(et) term only influences a simple harmonic oscillator to O(e) 
on time scales of O(-) - expresses a weak averaging effect, and can be proved by integration by parts. Of 
course, this argument can be modified to give information about perturbed oscillators on longer times scales 
of O(t-), a < 2, and many different generalizations are possible. 

A simple generalization, which is usful for the study of slow motion of solitons, can be obtained by 
integrating the identity 



where F, G are sufficiently regular functions of t, x but u = u(t) depends on t only and the inner product is 
L 2 (dx). This is often useful because in perturbation theory for solitons functions often arise with (d t +u- V)G 
small - see ([T3g|) . 

The following result, used in the proof of proposition 1241 is a more complicated version of this idea: 
Proposition 40 In the situation of lemma \23[ 



where a function f satisfies f = d/ dt ( o ( W ) J if there exists a C 1 function g = o I W j and f = -^g. 



Proof We work mostly with the potential Vi(</>) = —U(\<j>\) which determines AT: recall that V[((j>) = 
— /3(\<p\)<p, and (being slightly cavalier with notation) (|76|) can be rewritten 




((fit + u • V)F, G} L 2 = d t (F, G) L 2 - (F, (d t +u- V)G) L2 



(195) 



((d t +ih + u- V)v,Af(f u , e J u ,v)) L3 =d t (o (V)) +0(eW + e 3 W^ 





mU,e, U,V) = -Vi(/ W , e +V)+ Vi (/„,«) + Vi (/„)(«) 



Define 




(196) 




Exp[i@]f L 



(197) 




(198) 



(d t v + ihv,J\f(f u>e , f u , v)) 



L2 — 
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But, 



L 2 



ft I (Vi(/;, e + «*)-v 1 (/; ie )-v;(/: ie )K]--vr(/:, e )K] 2 



dx 



(200) 



L- 



Hence, 



(d t v*,V{(ti je + «*) - V[(f* >e ) - V{{fZ, e )[v*]) L2 = 

ft J (Vi(f u , a +v)- ViCL, e ) - V[(f u , e )[v] - (/ W) e)M 2 ) 

- ((ft + ifc) A,, e , V((/ W , e + «) - V((/ W , e ) - V'l{U, e )[v]) L2 - 



dx 



(dt + ih)f u , e ,^v[ 3) (U, e )[v} 2 



L' 2 



Now, 



V| 8) (/ W , e )[«] 2 )^ <C j \f u , e h\ (l + |/^ e | 3 ) \v\ 2 d X , 

by condition (fT3|l . Additionally, 

(i/i/^e,V(CL, e + «) - V((/ W , e ) - W„,e) M) L 2 
L 

(1 - fl) {ihf u , e , (V?(f u , e + SV) - Vr(/w,e)) M) L 2 , 
<C / \f u , e h\ (l + |/^ e | 3 ) (\v\ 2 + \v\ 5 )dx, 



(201) 

(202) 
(203) 



by condition (fT3)) . Therefore, by the exponential decay of f^.e and the fact that \f u ,eh\LP = 0(e) by the 
bounds of g2TH 



ifc/ U) «Vl(/ W)e + «) - Vi(f u , e ) - (/„,»)[«] - ^Vf 3) (/ W , e )[ 



= 0(eW). 



L 2 



Integration by parts and lemma [35] imply that 
((u • Vv, A/"(/ u , e , f u ,o,v)) = 0(e 2 W)+ 

(u ■ V/ UlOJ V[{U e + v)- V[(f u>e ) - V'{{U, e )[v] - \v?\f u ,*)[vA . 

Next notice that the quantity 

/(ft + u • V)/„, e , VU/„, e + v)- V[(f u , e ) ~ V{ (f u , e )[v] - Iv^if^vA 

can be estimated to be O(eW) in the same way as the bounds (|202[1 . (|203[) once we note that, for every 
P G [l,oo], 



||ft/ w , e + u-V/ w 



LP 



O(e), 



(204) 



(A - V (X)) ■ d x f u 

by (]102[) . The proof is now completed by noticing that Taylor's theorem and (|T3|) imply that the quantity 

/ (vi(/:, e + v*) - vixj - w;,>*] - ^'(.cjh 2 ) & 

is 0(W 3 / 2 ) + 0(W 3 ) = o(W), since W is small by assumption. 



□ 
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